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Abstract 

The Relativistic Continuum Hartree-Bogoliubov (RCHB) theory, which properly takes into account the 
pairing correlation and the coupling to (discretized) continuum via Bogoliubov transformation in a micro- 
scopic and self-consistent way, has been reviewed together with its new interpretation of the halo phenomena 
observed in light nuclei as the scattering of particle pairs into the continuum, the prediction of the exotic 
phenomena — giant halos in nuclei near neutron drip line, the reproduction of interaction cross sections and 
charge-changing cross sections in light exotic nuclei in combination with the Glauber theory, better restora- 
tion of pseudospin symmetry in exotic nuclei, predictions of exotic phenomena in hyper nuclei, and new 
magic numbers in superheavy nuclei, etc. Recent investigations on new effective interactions, the density 
dependence of the interaction strengthes, the RMF theory on the Woods-Saxon basis, the single particle 
resonant states, and the resonant BCS (rBCS) method for the pairing correlation, etc. are also presented in 
some details. 
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1 Introduction 

Currently nuclear physics is undergoing a renaissance as evidenced by the fact that worldwide there are lots 
of Radioactive Ion Beam (RIB) [1] facilities operating, being upgraded, under construction or planning to be 
constructed, e.g., the Cooling Storage Ring (CSR) project at HIRFL in China will be completed in 2005 [2], 
the RIB Factory at RIKEN in Japan will begin to operate at the end of 2006 [3], the FAIR project at GSI in 
Germany was approved in 2003 [4], and the construction of the Rare Isotope Accelerator (RIA) is considered 
to be of the highest priority of all physics disciplines in the US [5], etc. These new facilities together with 
developments in the detection techniques have changed the nuclear physics scenario. It has come into reality 
to produce and study the nuclei far away from the stability line — so called "EXOTIC NUCLEI" [6-12]. 
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New exciting discoveries have been made by exploring hitherto inaccessible regions in the nuclear chart [6- 
12]. For examples, a new phenomenon, halo — a state in which nucleons spread like a thin mist around the 
nucleus, was first discovered in ^^Li with RIB in 1985 [13] and later in many other light exotic nuclei as well; 
exotic nuclei exhibit other interesting phenomena such as the disappearance of traditional shell gaps and the 
occurrence of new ones, which result in new magic numbers [14]; etc. The exotic nuclei also play important 
roles in nuclear astrophysics as well, as their properties are crucial to understand stellar nucleosynthesis. 

The current nuclear models are mainly based on the knowledge obtained from the nuclei near the /?- 
stability line. For instances, the cornerstones for the edifice of modern nuclear physics include shell model [15] 
and collective model [16, 17], which are respectively based on the magic numbers (the stability of some nuclei 
compared with their neighbors) and the incompressibility of the nuclear matter. Therefore the change of 
magic numbers and the unprecedented low density nuclear matter in halo nuclei have shaken the foundation 
of nuclear physics. New innovative nuclear models are needed to describe the exotic nuclei characterized by 
weakly binding and low density. 

The relativistic mean field (RMF) theory [18] has received wide attention due to its successful description 
of lots of nuclear phenomena during the past years [19-23]. In the framework of the RMF theory, the nucleons 
interact via the exchanges of mesons and photons. The representations with large scalar and vector fields in 
nuclei, of order of a few hundred MeV, provide simpler and more efficient descriptions than nonrelativistic 
approaches that hide these scales. The dominant evidence is the spin-orbit splittings. Other evidence includes 
the density dependence of optical potential, the observation of approximate pseudospin symmetry, correlated 
two-pion exchange strength, QCD sum rules, and more [24]. The relativistic Bruekner-Hartree-Fock theory 
and RMF theory with density-dependent coupling constants extracted from it can reproduce well the nuclear 
saturation properties (the Coester line) in nuclear matter [25,26]. Furthermore the RMF theory can reproduce 
better the measurements of the isotopic shifts in the Pb-region [27], give more naturally the spin-orbit potential, 
the origin of the pseudospin symmetry [28, 29] as a relativistic symmetry [30-32] and spin symmetry in the 
anti-nucleon spectrum [33], and be more reliable for nuclei far away from the line of /3-stability, etc. Obviously, 
RMF is one of the best candidates for the description of exotic nuclei. 

In order to describe exotic nuclei, the pairing correlation and the coupling to continuum, which are extremely 
crucial for the description of drip line nuclei, must be taken into account properly [34]. The continuum effect 
is commonly taken into account in the Hartree-Fock-Bogoliubov (HFB) [35] or relativistic-Hartree-Bogoliubov 
(RHB) [21] approaches. In most of these calculations the continuum is replaced by a set of positive energy states 
determined by solving the HFB or RHB equations in coordinate space and with box boundary conditions [36- 
38]. Recently the HFB equations were also solved with exact boundary conditions for the continuum spectrum, 
both for a zero range [39] and a finite range pairing forces [40]. These resonant continuum HFB or RHB 
results with exact boundary conditions are generally close to those with box boundary conditions [39,41]. The 
extension of the RMF theory to take into account both bound states and (discretized) continuum via Bogoliubov 
transformation in a microscopic and self-consistent way, i.e., the Relativistic Continuum Hartree-Bogoliubov 
(RCHB) theory has been done in Refs. [36,42]. The RCHB theory was very successful in describing the ground 
state properties of nuclei both near and far from the /3-stability line. By using a density-dependent zero 
range interaction, the halo in ^^Li has been successfully reproduced in this self-consistent picture. Remarkable 
successes of the RCHB theory include the new interpretation of the halo in ^^Li [36], the prediction of the 
exotic phenomena as giant halos in Zr {A > 122) [43] and Ca {A > 60) [44] isotopes, the reproduction of 
interaction cross section and charge-changing cross sections in light exotic nuclei in combination with the 
Glauber theory [45,46], better restoration of pseudospin symmetry in exotic nuclei [31,32], and predictions of 
exotic phenomena in hyper nuclei [47] and new magic numbers in superheavy nuclei [48], etc. 

The main purpose of the present manuscript is to review the RCHB theory and its applications for exotic 
and superheavy nuclei in spherical case. Some other related topics will also be covered briefiy. In Section [21 
the formalism, the numerical solutions and the effective interactions for the RMF theory are presented. In 
Section 13 the pairing correlations and the approaches for pairing are sketched. Discussion on the continuum 
states follows where the resonant BCS method and several methods to obtain single particle resonant states 
are briefiy reviewed. Then the detailed formalism of the RCHB theory together with the discussion on the 
effective pairing interactions are given. In Section [IJ the applications of the RCHB theory to the properties 
of exotic nuclei are presented, such as the binding energies, particle separation energies, the radii and cross 
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sections, the single particle levels, shell structure, the restoration of the pseudo-spin symmetry, the halo and 
giant halo, and halos in hyper nuclei, etc. In Section 13 the predictions of new magic numbers in superheavy 
nuclei are presented. Finally a brief summary and perspectives are given in Section El 



2 Relativistic mean field theory 

In this section, we present the formalism, the numerical solutions and the effective interactions for the RMF 
theory as well as its application for nuclear matter. In the first subsection, the effective Lagrangian density and 
equations of motion for the nucleon and the mesons are given. The numerical solutions of the Dirac equation 
for finite nuclei are discussed in the second subsection. Subsequently, the effective interactions with nonlinear 
self-coupling meson fields and density dependent meson-nucleon couplings, and their influences on properties 
of nuclear matter are discussed. 



2.1 The general formalism 

The RMF theory describes the nucleus as a system of Dirac nucleons which interact in a relativistic covariant 
manner via meson fields. The meson fields are treated as classical fields. In the simplest RMF version, i.e., 
the a-uj model [18], the mesons do not interact among themselves, which leads to too large incompressibility 
in nuclear matter. Therefore a nonlinear self-coupling of the cj-field was proposed [49]. In order to repro- 
duce the density dependence of the vector and scalar potentials of the Dirac-Brueckner calculations [26], the 
nonlinear self-coupling of the w-meson was found to be necessary [50]. Recently the nonlinear self-coupling of 
the isovector p-meson was also introduced to improve the density-dependence of the isospin-dependent part of 
the potentials [51]. Within the present scheme, the isoscalar-scalar u-meson provides the mid- and long-range 
attractive part of the nuclear interaction whereas the short-range repulsive part is provided by the isoscalar- 
vector w-meson. The photon field A^^{x) accounts for the Coulomb interaction while the isospin dependence 
of the nuclear force is described by the isovector-vector p-meson. The vr-meson field is not included because it 
does not contribute at the Hartree level. In principle, other mesons apart from a, uj, and p may also contribute 
to the nuclear interaction as well, e.g., the isovector scalar (5-meson, which was suggested in Ref. [52] and also 
on the basis of a relativistic Brueckner theory in Refs [53-56], has been used in some structure calculations for 
exotic nuclei [55], nuclear matter [57] and stellar matter [58]. However, as the RMF theory is only an effective 
theory therefore it is expected that the contributions from other mesons can be effectively taken into account 
by adjusting the model parameters to the properties of nuclear matter and finite nuclei. With various versions 
of the nonlinear self-couplings of meson fields, the RMF theory has been used to describe lots of nuclear phe- 
nomena during the past years with great successes [19-21]. In order to avoid the problem of instability for 
the nonlinear interactions at high densities, the RMF theories with density dependence in the couplings are 
developed as well [26,51,59-62]. 

The Lagrangian density of the RMF theory can be written as 

l-rs" 



+ ^^d^'ad^a - U^{a) - Ul^^'^^u + - \r'''' ■ R^u + Up{p^) (1) 

where M and mi{gi) (i = a,uj,p) in the following are the masses (coupling constants) of the nucleon and the 
mesons respectively and 

nt^'^ = d^'u;" -d^Lof", (2a) 
= dt'^-d"^, (2b) 
pt^u ^ Qt^^i^ _Qyj^t^ (2c) 

are the field tensors of the vector mesons and the electromagnetic field. We adopt the arrows to indicate vectors 
in isospin space and bold types for the space vectors. Greek indices /i and v run over 0, 1, 2, 3 or t, rc, y, z, 
while Roman indices i, j, etc. denote the spatial components. 
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The nonlinear self-coupling terms Ua{cr), U^{uj^)^ and Up{p^) for the fi-meson, w-meson, and p-meson in 
the Lagrangian density respectively have the following forms: 



1 O 2 Q 1 A 

+ 3520-^ + -53(7 , 



(3a) 
(3b) 
(3c) 



2 P'- '-^ ■ 4 

From the Lagrangian density the Hamiltonian operator can be obtained by the general Legendre 
transformation 

dcl)i{x) 



H 



(fx n 



(fx 



^T:i{x)- 



dt 



C(x) 



where the conjugate momenta of the field operators 

TTi{x) = 



ilj^a^uJu-^Pu^-^u) are defined as 



dC 



d [d(t>,/dt] ■ 

Then the Hamiltonian density of the system can be easily obtained as 



(4) 



(5) 



-• - - - -1 — Ts 



(6) 



-n ■ + -n ■ Uii. 
4 4 



Upipp 



The relativistic mean field theory is formulated on the basis of the effective Lagrangian with the mean 
field approximation, i.e. the meson fields are treated as classical c-numbers. With this approximation, all the 
quantum fluctuations of the meson fields are removed and the nucleons are described as independent particles 
moving in the effective meson and photon fields. Therefore the nucleon field operator can be expanded on a 
complete set of single-particle states as, 

V'(2;) = ^V'a(x)Ca, (7) 
a 

where Ca is the annihilation operator for a nucleon in the state a of the Dirac fields and ipa is the corresponding 
single-particle spinor. The operator Ca and its conjugate Ca satisfy the anticommutation rules 



{ca, 4} = Sab and {ca,Cb} = {cjj, cl} = 0. 



(8) 



Confined to the single-particle states i with positive energies, i.e. the no-sea approximation, the ground state 
of the nucleus can be constructed as, 



|$) = JJ4|0) with ($1$) = 1, 



(9) 



i=l 



where |0) is physical vacuum. 

With the ground state © and the mean field approximation, the energy functional, i.e. the expectation 
value for the Hamiltonian © is obtained as 

(fxTr 



+ / d-^x 



/? (^7 • p + M + + (^..w'^Tm + 9pr ■ P^,l^ + -e(l - n)Ap-i>'j p 



(10) 
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where (j) = a,uju, p^, Ai, and the density matrix p is defined as 

Pij = (^IV'JcJViCil^) = vJV'i'Ju- (11) 

For system with time reversal symmetry, the space-like components of the vector fields vanish. Furthermore 
one can assume that in all nuclear applications the nucleon single-particle states do not mix isospin, i.e. the 
single-particle states are eigenstates of t^, therefore only the third component of p^, survives. Stationarity 
implies that the nucleon single-particle wave function can be written as 

Mx) = e-'^'Vi(x). (12) 

where is the single-particle energy. Accordingly the density matrix is reduced as, 

= '0](x')'0j(x)(5ij. (13) 

Altogether there remain only the meson fields a, loq, /9q, Aq which are time independent. 

The equations of motion for the nucleon and the mesons can be obtained by requiring that the energy 
functional (|lflj) be stationary with respect to the variations of p and (p. More explicitly, the stationary condition 
reads 

S{EnMF{p,(l))-Tv{£p)) = 0, (14) 

where f is a diagonal matrix, whose diagonal elements are the single particle energies ei{i = 1, • • • , N) intro- 
duced in Eq. (|12|) . and N is the number of eigenstates. 

Using the variation dp with respect to tpi, the stationary condition 1)14(1 leads to the Dirac equation 

(a • p + (3{M + S) + V) V'i(x) = eiVi(x) (15) 

for the nucleon and the Klein-Gordon equations for sigma, omega, rho, and the photon: 

-VV + [/;(a) = -g^p,, (16a) 

-V'^ujo + U^{uJo) = QuiPv, (16b) 

-V^pl + U'p{pl) = Qpp^, (16c) 

-V^Ao = epc. (16d) 
The scalar potential S and vector potential V in equation ()15|) are respectively: 

5'(x) = gaO-(x), (17a) 

F(x) = 5<^cuo(x)+5pT3p|](x) + ^e(l-r3)^o(x). (17b) 

While the scalar density {ps), the baryon density {pv), the isovector density (/O3), and the charge density (pc) 
in the Klein-Gordon equations (|16aH16d|) are respectively, 

Ps = Tr[/3p], (18a) 
Pv = TV[p], (18b) 
P3 = TV [rap], (18c) 

Pc = ^Tr[(l-T3)p]. (18d) 



The total energy of the system can be obtained from the energy functional (|1U() as, 

E= j Sx |Tr [(a • p + /3M) p\ + ^Tr g„l3ap + g^i^op + gpnplp + ^e(l - T3)^o 



(19) 



In the density dependent RMF approach, where the nonlinear self-couplings for the a, lo, and p mesons in the 
Lagrangian density are respectively replaced by the density dependence of the coupling constants g^ip), duiip), 
and gp{p), an additional term, i.e. the rearrangement term, will appear in the Dirac equation ((T3|l [26,51,59-62]. 
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2.2 Numerical algorithm for spherical nuclei 

The harmonic oscillator basis has served as a very useful tool in nuclear structure study. Normally, the equations 
of motion for nucleons moving in a mean field are solved by expanding them on the harmonic oscillator (HO) 
basis [15-17,35,63]. However, for exotic nuclei with large spacial extension, e.g., halo nuclei, it is not justified 
to work in the conventional harmonic oscillator basis due to its localization [36,64-66]. Instead, one can choose 
to work either in the coordinate space, or improve the asymptotic behavior of the HO wave function, or adopt 
other basises which have a correct asymptotic behavior, for example, the Woods-Saxon basis. 

In this subsection, we will focus on the numerical solution of the RMF for spherical nuclei. Due to the 
special spacial symmetry, both the Dirac equation for the nucleon and the Klein-Gordon equations for the 
mesons and the photon become radially dependent only, thus facilitating much the solution of the coupled 
equations. The formalism for the spherical relativistic Hartree (SRH) theory will be briefly presented. We 
then review the SRH theory in different basis, including Finite Element Method (FEM) in coordinate space, 
transformed harmonic oscillator basis, and the Woods-Saxon basis. The application of the SRH theory to 
doubly magic nuclei follows. 

2.2.1 Spherical relativistic Hartree theory 

Starting from the Eqs. (jlSj) and (|16aH16d)) given in the previous subsection, one derives the coupled radial 
equations for spherical nuclei, i.e., the radial Dirac equation and radial Klein-Gordon equations. 

For spherical nuclei, the Dirac spinor which is the expansion coefficient V'a(r) (a = {a,K, m}) in Eq. (|12j) 
(the coordinate from x has been changed to r to reflect the spherical symmetry) is characterized by the angular 
momentum quantum numbers K{l,j), m, the parity, the isospin t = ±1/2 ("+" for neutrons and "— " for 
protons) and the radial quantum number a and has the form 



with G'^{r)/r and F^{r)/r the radial wave functions for the upper and lower components and Yj^{9,(j)) the 
spinor spherical harmonics [67]. Substituting Eq. H2()j) into the Dirac equation (|15|) . one can deduce the radial 
Dirac equations as 




(20) 




(21b) 




with the scalar and vector potentials 




(22a) 
(22b) 



The meson field equations (|16aH16cH) simply become radial Laplace equations of the form 




(23) 



where are the meson masses for cf) 



cr,uj,p and zero for the photon. The source terms are 




for the cr— field, 
for the w— field, 
for the p— field, 



for the Coulomb field, 



(24) 
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with 



1 ^ 



47]- 7-2 

2=1 

1 ^ 

/^-W = i;^E(l^^(^)l' + I^^Wl')' (25b) 

2 = 1 
1 ^ 

^3(0 = 4;^E^3(|G.(r)|2 + |F,(r)|2), (25c) 

2 = 1 

'^^(^) = 4;^E^(l-^3)(|G'2(r)p + |F,(r)|2). (25d) 

2 = 1 

The procedures of solving these coupled equations are as the following: a) with a set of estimated meson 
and photon fields, the scalar and vector potentials in Eas. ()22() are calculated and the radial Dirac equation 
is solved; b) the so obtained nucleon wave functions are used to give the source term of each radial Laplace 
equation for mesons and the photon; and c) the new meson and photon fields obtained from solving these 
Laplace equations will be used to replace the fields in step a). This procedure is iterated until a demanded 
accuracy is achieved. 

With the above procedures, there are two methods to solve the coupled equations (PIIE^ . One is in 
coordinate space with the shooting method [68], i.e. spherical relativistic Hartree theory in r space (SRHR), 
and the finite element method [42,69-71] (SRHFEM). Another is in configuration space, e.g., the harmonic 
oscillator basis [72] (SRHHO), the transformed harmonic oscillator basis [64] (SRHTHO) and the Woods-Saxon 
basis [73] (SRHWS). The readers are referred to Ref. [68] and Ref. [72] for SRHR and SRHHO, respectively. 
In the following the numerical techniques for the solution of the Dirac equation in the finite element method, 
the transformed harmonic oscillator basis and Woods-Saxon basis will be briefiy introduced. 

2.2.2 Finite element method 

A convenient procedure for the coordinate space discretization of the Dirac equations 1)21(1 is provided by Finite 
Element Method (FEM) [42,69-71]. Similar as in the r space [68], the Dirac equations (|21|) are solved in a box 
with a box size R and proper boundary condition. The radial wave functions and are discretized at 
A'^ + 1 points, i.e., ri = 0, r2 = Ar, ... , rjv+i = N ■ Ar, where Ar = R/N . The wave functions thus obtained 
become tabulated values: G%{rQ),G%{ri), ... , G^{rN+i) and F^(ro),F^(ri), ... , F^{rN+i)- 

The idea of FEM is as following: if the Ar = R/N is small enough, or equivalently N is large enough, the 
wave function G'^{r) in the element r G [rj,rj+i] can be well approximated by G'^{ri) and GJJ(rj+i) together 
with some simple analytic function, e.g., the linear, quadratic, cubic, or 4th order shape function [69]. 

Taking the linear shape function, 

' - ^' (26) 
P, 

with p = — - — - and r S [ri,rij^i] as an example, the wave function GJ^(r) can be well approximated by. 




Ar 



GS(r) = GS(r,)r?a(r) + GS(r,+i)%(r), r G [r„r,+i]. (27) 



Similarly the expression for F^{r) can also be obtained. For r G [r^, rj+i], the Dirac equations (|21|) can now be 
written as: 



M + S{r) + V{r)-e^ + - ^ 

ar r 

^ + - -M-S{r) + V{r)-e^ , 

ar r I 



G'S(r2)r/a(r) + G«(ri+i)j?fc(r) 
^a(?'2)??a(r)+F^(ri+i)j76(r) 



(28) 
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Multiplying (G^(rj)r/a(r)+G;j(rj+i)77fe(r) F^(rj)77a(r)+F^(rj+i)r?fe(r)) from left and integrating with / r dr, 

J Vi 

algebra equations for GJJ(ri), GJJ(rj+i) and F^{ri), F^(rj+i) are obtained. Repeating the same procedure for 
all the elements, the Dirac equations (j2^ are transformed into a generalized eigenvalue problem. This 
generalized eigenvalue problem can be solved by the standard diagonaliztion algorithms and all the energies 
and wave functions discretized at r = ri, r2, tat+i can be obtained. It has been shown that FEM can 
provide very accurate solutions for the relativistic eigenvalue problem in the self-consistent mean-field approx- 
imation [42,69-71]. 

Comparing with the shooting method in the r space [68], the FEM has the advantage that one can get 
all the energy Cq^ and wave functions and by a single diagonalization and it is straight forward to be 
generalized to the cases with nonlocal interaction in the pairing channel. However, in order to obtain the same 
accuracy as the shooting method, a huge matrix has to be constructed and it becomes more time consuming. 
One can also replace the linear shape function ria{r) and ryf,(r) in Eq. (|26|) by the quadratic, cubic, or 4th order 
shape function, the procedures are the same and in principle the same accuracy can be achieved by the linear 
shape function with larger A^. 



2.2.3 Transformed harmonic oscillator basis 

In order to modify the asymptotic behavior of the harmonic oscillator wave function at large r, the local- 
scaling transformation method [74-76] has been introduced to construct the so called local-scaling transformed 
harmonic oscillator basis (THO) in Refs. [64,65]. 

A local-scaling point coordinate transformation (LST) is defined as 

r' = f(r) ^ f-/(r). (29) 

The transformed radius vector has the same direction f = r/|r|, while its magnitude r' = f{r) depends on the 
scalar LST function, /(r) is assumed to be an increasing function of r, and /(0)=0. The corresponding LST 
wave function can be expressed as 



^/(ri,r2,...,r^) 



n 



/^(r,)8/(r,- 



^,2 

.4 = 1 



'l'(/(ri),/(r2),...,/(rA)), (30) 



where ^(ri, r2, r^^) is an A-particle wave function normalized to unity. The local one-body density corre- 
sponding to an A-body wave function \I' is 

pir)=A j |^'(r,r2,...,r^)|2dr2...drA. (31) 

There exists a simple relation between the local density p/(r) associated with the LST function and the 
density p{r) which corresponds to the prototypical model function 'J': 

f{v) dfjr) 

Pfi"") = —2 ^/'(/W)- (32) 

When the form of the density p/(r) is known, Eq. ()32() becomes a first order nonlinear differential equation for 
the LST function / and can be solved easily. For a system with spherical symmetry, Pf,p, and / depend only 
on r=|r|, and Eq. H32|) can be reduced to a nonlinear algebraic equation 

r rfir) 

Pf{u)u'^du = I p{u)u^du. (33) 
Jo 

The solution can be found subject to the boundary condition /(O) = 0. 

For shell-model or mean-field applications, one has to consider the case when the model many-body wave 
function is a Slater determinant 

^'(ri, r2, r^) = det |^i(r,)|. (34) 
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Ne coordinate calculations 
THO calculations 
HO calculations 




Figure 1: Self-consistent RHB proton and neutron densities for the ground-states of ^^Ne and ^^Ne. Density 
profiles calculated in THO are compared with FEM. Taken from Stoitsov et al. [64]. 

The single-particle functions 0i(r) form a complete set. The LST wave function is defined by the transformation 
pn]) . and is written as a product state 



*/(ri,r2, ■■■,rA) 



det |(/)i(rj 



of the transformed basis state 



f\r) dfjr) 

J.2 



1/2 



urn)- 



(35) 



(36) 



In Refs. [64,65], the transformed harmonic oscillator basis (THO) has been derived by a local scaling-point 
transformation of the spherical harmonic-oscillator radial wave functions. The unitary scaling transformation 
produces a basis with improved asymptotic properties. The THO basis is employed in the solution of the 
relativistic Hartree-Bogoliubov (RHB) equations in configurational space [64]. As shown in Fig.^ an expansion 
of nucleon spinors and mean-field potentials in the THO basis reproduces the asymptotic properties of neutron 
densities calculated by FEM in the coordinate space. 



2.2.4 Woods-Saxon basis 

Woods-Saxon basis from the Schrodinger equation (SWS basis) For the Schrodinger equation with 
a spherical Woods-Saxon potential 



Vo 



r <Rn 



I -)- g{r-Ro)/cio ' 



(37) 



where i?max is introduced for practical reasons to define the box boundary. The eigenfunction can be written 
as (pnimii^) = Rni{r)Yimi{d , (p) and its radial Schrodinger equation is derived as 



2M V r'^ dr dr 



Rnl{r) = EnlRnl{r). 



(38) 



Equation (|38|) is solved on a discretized radial mesh with a mesh size Ar. i?max (Ar) is chosen large (small) 
enough to make sure that the final results do not depend on it. The radial wave functions thus obtained form 
a complete basis, 

{i2„z(r);n = 0,l,--- ;/ = 0,l,--- ,n}, (39) 
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in terms of which the radial part of the upper and the lower components of the Dirac spinor in Eq. (|21() are 
expanded respectively as 



n=0 

^max 



n=0 



The radial Dirac equation (PT|) is transformed into the WS basis as 



^mn >~'mn 

V- - 

--mn ^mn 



9an 
fan 



9an 
fan 



where the matrix elements are calculated as follows 



•^mn 

E - 

C-mn 



r'drR^l{r) {V{r) + S{r) + M) R^i{r), 



r drRmi{r) ( + — 



Rn 



r 



r^drR^^l{r) 



d + 1 



dr 



Rnlir), 



r^drR^jir) {V{r) - S{r) - M) R^j{r). 



(40) 



(41) 



(42) 



In practical calculations, an energy cutoff Ecut (relative to the nucleon mass M) is used to determine the 
cutoff of the radial quantum number nmax [73]. 



Woods-Saxon basis from the Dirac equation (DWS basis) The radial Dirac equation (|21() may be 

solved in r space [68] with Woods-Saxon-like potentials for Vo{r) ± So{r) [77] within a spherical box of the size 
-Rmax) together with the spherical spinor which gives a complete WS basis 



i [e° ib'^ (r s t)] ■ e° < O'^ 



(43) 



with n = 0, 1, • • • , K = ±1, ±2, • • • , and m = —Jk, ■ ■ ■ , Jk- "^nKmi^^ takes the form of Eq. l|20|) . In such 
cases, states both in the Fermi sea and in the Dirac sea should be included in the basis for completeness. The 
nucleon wave function (|2flj) can be expanded in terms of this set of basis as 



V'oKm(r, S, t) ^ ^ ^an'4'nK,mi^' ^^^)^ 



(44) 



n=0 



where rir. 



'^max + '^max + ^ ^nd the Summation runs over positive energy levels in the Fermi sea for 



^max '"max ' '"max 

< n < and over negative energy levels in the Dirac sea for -|- 1 < re < remax- The negative energy 
states are obtained with the same method as the positive energy ones. In this WS basis, the Dirac equation 
(PT|) turns out to be 



(45) 



n=0 



with 



HLn = (V'^(r)|[Ay(r) + /5A5(r)]|^0(r)> 



+ 



Rn 



Rn 



drG'^ir) [AV{r) + AS{r)] G°(r) 
drF^ir) [AVir)-AS{r)]F^{r), 



(46) 
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Table 1: Effects of negative energy levels on bulk properties in SRHDWS for O. The effective interaction 
for the Lagrangian is NLSH, .Rmax = 20 fm, Ar = 0.1 fm, and N^^^ = 25. For the initial Woods-Saxon like 
potentials, parameters in Ref. [77] are used except for Vq which is specified in the table. The left value in each 
entry gives the result without negative energy levels included and the right one that with iV~3^ = 25. Energy 
is in MeV and radius in fm. Taken form Ref. [73]. 



Vo 


54 MeV 


72 MeV 


90 MeV 


N~ 

'max 





25 





25 





25 


E/A 


8.547 


8.013 


8.117 


8.015 


8.427 


8.012 




2.385 


2.568 


2.531 


2.567 


2.610 


2.567 



where AV{r) = V{r) — Vo(r), AS{r) = 5(r) — So{r) and the angular, spin, and isospin quantum numbers are 
omitted for brevity. 

In the expansion (|44|) of the nucleon wave function in the SRHDWS theory, one has to take into account 
not only the states in the Fermi sea but also those in the Dirac sea because these states form a complete basis 
together. The contribution from negative energy states for ^^O is given in Tabled It is found that, without 
including the negative energy levels, the calculated results will depend on the potentials for the basis. 

The contribution of negative energy states in the Dirac sea to the wave function can be calculated by 
X^n I'-n P the expansion H44|) and it is around 10~^~~^ (Note that the nucleon wave function is normalized 
to one). However, such a small component from negative energy states in the wave functions contributes to 
the physical observables such as E/A and rrms by magnitudes of 1~10 % as can be seen from Tabled 

Comparisons between r space, harmonic oscillator basis and Woods-Saxon basis It is found that 
for stable nuclei, the SRHR, SRHSWS, SRHDWS, and SRHHO approaches are all valid and results from them 
are in excellent agreement with each other [73]. However, for unstable nuclei near the neutron drip line, these 
methods differ from each other to some extent. 

In Fig. 121 the neutron density distribution of ^^Ca from different SRH approaches are compared. With 
the same box size, the density distribution from SRHR are almost identical with those from SRHWS, which 
indicates the equivalence between SRHWS and SRHR. For brevity, only pn{r) from SRHR with i?max = 35 fm 
is displayed which covers the curve corresponding to Pn{r) from SRHWS with i?max = 35 fm in Fig. [3 On the 
other hand, Pnir) from SRHHO even with A^n^ax = 43 fails to reproduce the result of SRHR due to the well 
known localization property of HO potential [64] . 

These results indicate that even the long tail (or halo) behavior in neutron density distribution for nuclei 
near the drip line can be reproduced quite well by the expansion method in the Woods-Saxon basis in the 
Schrodinger framework (SRHSWS) and that in the Dirac framework (SRHDWS). This can be seen that the 
neutron density distribution obtained in SRHR is reproduced well in both SRHSWS and SRHDWS when 
enough box size is taken. 

2.3 Effective interactions 

In the Lagrangian density there are meson masses m„,mi^,mp and meson-nucleon coupling constants 
g<j,guj,gp together with the nonlinear self-couplings of the meson fields left to be determined. They are the 
nucleon-nucleon interactions in the RMF theory and in principle should be determined either by more funda- 
mental theories or by experiments. However, as the relativistic mean field theory is formulated on the basis of 
the above effective Lagrangian in connection with the mean-field and the no-sea approximations, it is difficult 
to determine these interactions microscopically. Instead, the masses and coupling strengthes of the mesons and 
the nonlinear self-couplings of the meson fields are determined by reproducing the properties of nuclear matter 
and a few doubly magic nuclei. Namely, they are effective interactions in the similar sense as their conventional 
counterparts. The effective interactions in the RMF theory could be determined by minimizing the least square 
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Figure 2: Comparison of density distributions for ^^Ca from SRHR, SRHWS and SRHHO with NLSH. Ar = 
0.1 fm for SRHR and SRHWS. ^cut = 75 MeV and fi^ax = 20 (thick dot-dashed curve), 25 (thick long-dashed 
curve), 30 (thick dashed curve) and 35 fm (thick dotted curve) for SRHWS. These sets of cutoff's correspond 
to cutoff's in principal quantum number A'max = 25 (thin dot-dashed curve), 31 (thin long-dashed curve), 37 
(thin dashed curve) and 43 (thin dotted curve) which are used in SRHHO calculations. Taken form Ref. [73]. 

error as follows, 

(47) 

where a is the ensemble of the meson masses rriu , , nip and the meson- nucleon coupling constants ga,9uj, Qp 
together with the nonlinear self-couplings of the meson fields to be fitted, and y^^^ and fjj are the experimental 
observables and corresponding weights. In general, the masses and charge radii of spherical nuclei near the 
/3-stability line are adopted as observables in the least-square fitting procedure. 

Among the existing effective interactions for the RMF theory, the most frequently used are NLl [78], 
NLSH [79], TMl [50] and NL3 [80] with nonlinear self-couplings of mesons. Along the /3-stability line NLl 
gives excellent results for binding energies and charge radii, in addition it provides an excellent description 
of the superdeformed bands [81,82]. However, when moving away from the stability line the results are less 
satisfactory. This can be partly attributed to the large asymmetry energy J ~ 44 MeV. Moreover, the neutron 
skin thicknesses calculated with NLl show systematic deviations from the experimental data [83]. 

In NLSH this problem was treated in a better way and the improved isovector properties have been obtained 
with an asymmetry energy of J ~ 36 MeV. Furthermore, NLSH seems to describe the deformation properties 
in a better way than NLl does. However, NLSH produces a slight over-binding along the line of /3-stability 
and in addition it fails to reproduce successfully the superdeformed minima in Hg-isotopes in constrained 
calculations for the energy landscape. A remarkable difference between the two effective interactions are the 
quite different values predicted for the nuclear matter incompressibility [84], i.e., K = 212 MeV for NLl while 
K = 355 MeV for NLSH [85,86]. As an improvement, NL3 and TMl provide reasonable compression modulus 
(i^NLS = 271.7 MeV,KTMi = 281.16 MeV) and asymmetry energy (Jnls = 37.42 MeV, Jtmi = 36.89 MeV) 
but fairly small baryonic saturation density (/Otmi = 0.145). In order to improve the description of these 
quantities, two nonlinear self-coupling effective interactions, PKl with nonlinear a- and w-meson self-couplings 
and PKIR with nonlinear cr-, uj- and p-meson self-couplings were developed [51] (see Table [21). 

In order to reproduce better the experimental quantities such as the binding energies and nuclear radii, etc., 
an additional correction should be added in calculating the energy, i.e., the center-of-mass correction. Conven- 




X^(a) = E 



Vi -y{xi]a.) 
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tionally a phenomenological center-of-mass correction as ——AlA is used. Microscopically the correction 
can be calculated by the projection-after- variation in first-order approximation [87]: 

i^=: = -2^(p?...>, (48) 

where the center-of-mass momentum Pc.m. = YlfPi the expectation value of its square (Pc.m.) reads 

(Pc.m. ) = Yl "^"^"^ ~ Y ^a^bPab 'Plb + Y ^'^^^^bUbPab ' Pab (49) 
a a,b a,b 

with occupation probabilities v'^ and u^^ = I — accounting for the pairing effects, where a, b denote the BCS 
states (see the following section). 

Here it should be mentioned that the prescription H48|) is based on non-relativistic considerations. It does 
not preserve Lorentz invariance. Furthermore, it also breaks the complete self-consistence of the variational 
scheme. However, as it is not included in the self-consistent procedure and only presents an additional correction 
term to the binding energy, one can be satisfied with it for the moment. Compared with the binding energy, 
this center-of-mass correction is sizable in light nuclei (about 9% in ^^O) but much less important in medium 
and heavy nuclei (about 0.4% in ^'^^Pb) as seen in Fig. |31 




16 32 48 64 80 96 112 128 144 160 176 192 208 

Mass Number 

Figure 3: The microscopic center-of-mass correction in comparison with two phenomenological cases. Taken 
from Ref. [51]. 

For the nuclear radii, the effects from the center-of-mass motion can also be taken into account as follows. 
Because of its fairly small effects, a rather rough correction is adopted for protons. 



6Rl 



Z 



2 

■cm. 



(50) 



where the center-of-mass coordinate Rr 



iMEfr*. Then one obtains, 



6Ri 



(51) 



where Rp and Rm denote the proton and matter radii. Here we only consider the direct-term contributions to 
Eq. (|5fl|) to conform to the spirit of the RMF theory. For the neutron radii, one can follow the same procedure 
as that for protons. The charge radius is obtained from the proton radius combining with the proton and 
neutron size and the center-of-mass correction (|51j) included in R^ [50] 



Rl^ = Rl + (0.862 fm)2 - (0.336 fmfN/Z. 



(52) 
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Table 2: The nonlinear effective interactions PKl, PKIR and density-dependent effective interaction PKDD 
[51] in comparison with TMl [50], NL3 [80], TW99 [61], and DD-MEl [62], 





PKl 


PKIR 


PKDD 


TMl 


NL3 


TW99 


DD-MEl 




939.5731 


939.5731 


939.5731 


938 


939 


939 


938.5000 


Mp 


938.2796 


938.2796 


938.2796 


938 


939 


939 


938.5000 




514.0891 


514.0873 


555.5112 


511.198 


508.1941 


550 


549.5255 




784.254 


784.2223 


783 


783 


782.501 


783 


783.0000 


rup 


763 


763 


763 


770 


763 


763 


763.0000 


9a 


10.3222 


10.3219 


10.7385 


10.0289 


10.2169 


10.7285 


10.4434 


9u, 


13.0131 


13.0134 


13.1476 


12.6139 


12.8675 


13.2902 


12.8939 


9p 


4.5297 


4.55 


4.2998 


4.6322 


4.4744 


3.661 


3.8053 


92 


-8.1688 


-8.1562 





-7.2325 


-10.4307 





0.0000 


93 


-9.9976 


-10.1984 





0.6183 


-28.8851 





0.0000 


C3 


55.636 


54.4459 





71.3075 








0.0000 


d3 





350 


















Tabic 3: Density-dependent parameters of PKDD [51] for meson-nucleon couplings in comparison with 
TW 99 [61] and DD-MEl [62]. 

ttcr C(j d(y ciijj bu di_j dp 

PKDD 1.327423 0.435126 0.691666 0.694210 1.342170 0.371167 0.611397 0.738376 0.183305 

TW99 1.365469 0.226061 0.409704 0.901995 1.402488 0.172577 0.344293 0.983955 0.515000 

DD-MEl 1.3854 0.9781 1.5342 0.4661 1.3879 0.8525 1.3566 0.4957 0.5008 
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With the microscopic center-of-mass motion, a multi-parameter fitting can be performed using the Levenberg- 
Marquardt method [88]. In Ref. [51], the masses of ^^O, "^^Ca, ^'^Ca, ^^Ni, ^^Ni, ^^Zr, ^^^Sn, ^^"^Sn, ^^''Pb and 
^'^^Pb and the bulk quantities of nuclear matter are chosen as observables to determine the effective interactions. 
The radii are excluded because the proper values of the compression modulus K and the baryonic saturation 
density po are sufficient to give a good description of the radii. For a fixed value of the compression modulus 
K, a large baryonic saturation density po will give a small charge radius and vice versa. Therefore a proper 
description of both masses and radii of finite nuclei could be obtained by carefully adjusting the values of these 
two quantities K and po- To give a fairly precise description of the masses, the center-of-mass correction is also 
essential for both light and heavy nuclei. As it can be seen in Fig. |31 the deviation between the microscopic 
and phenomenological results is considerably large not only for the light nuclei but also for the heavy ones. 
And there exist very remarkable shell effects in the microscopic results which are impossible to obtain with the 
phenomenological methods. The microscopic center-of-mass correction [87], therefore, is chosen to deal with 
the center-of-mass motion. 

Because the contribution to the nuclear masses from the nonlinear p-meson term is found to be fairly 
small, the effective interaction PKIR is obtained by fixing the nonlinear self-coupling constant ^3 to 350.0 and 
adjusting other parameters. 

In the RMF theory with density-dependent meson-nucleon couplings, the density-dependence of the coupling 
constants and can be parameterized as, 

9iiPv) = gi{Psai)fi{x) for i = a,uj, (53) 

where 

l + hi{x + diY 
^'^^^ = "- l + c.(x + d.)^ 

is a function oi x = p^/psati and psat denotes the baryonic saturation density of nuclear matter. For the p 
meson, an exponential dependence is utilized as 

9p = 5'p(Psat) exp[-ap(x - 1)]. (55) 

For the functions fi{x), one has five constraint conditions /i(l) = 1,/"(1) = /"(I) and //'(O) = 0. Then 
8 parameters related to density dependence for u-N and lo-N couplings are reduced to 3 free parameters. In 
general, the masses of the nucleons and the p- meson are fixed and the nonlinear self-coupling constants g2, 53, C3 
and ds are set to zero. With 4 free parameters for density dependence, totally there are 8~9 parameters left 
free in the Lagrangian density for the density-dependent meson-nucleon coupling RMF theory. A density- 
dependent meson-nucleon coupling effective interaction PKDD has also been obtained in Ref. [51] (see Tables |2l 
andlSJ. 

Tables 121 and El tabulate the new effective interactions PKl, PKIR and PKDD [51] in comparison with the 
old ones TMl [50], NL3 [80], TW99 [61] and DD-MEl [62]. The newly obtained ones reproduce better the 
experimental masses [89]. PKl, PKIR and PKDD also describe the charge radii very well, especially for those of 
the Pb isotopes. More comprehensive comparisons between Hartree-Fock-Bogoliubov, extended Thomas-Fermi 
model with Strutinski integral, RMF, and macroscopic-microscopic approaches with different forces have been 
performed for the description of nuclear masses and charge radii of spherical even-even nuclei (116 nuclides), 
from light (A=16) to heavy (A=220) ones in Ref. [90]. 

Tabled lists the nuclear matter quantities calculated with the newly obtained effective interactions PKl, 
PKIR and PKDD, in comparison with those from the other interactions. All the new effective interactions 
give a proper value for the compression modulus K. 



2.4 Density and isospin dependence of effective interactions 

There are so far quite a number of effective interactions, PKl, PKIR, PKDD [51] together with NLl, NL2 [91], 
NL3 [80], NLSH [79], TMl, TM2 [50], GL-97 [92] and the density-dependent effective interactions TW-99 [61], 
DD-MEl [62], etc.. It is very interesting to investigate the density and isospin dependence of the interaction 
strengthes of various effective interactions in the RMF theory and study their effects on nuclear matter [93]. 
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Table 4: Nuclear matter properties calculated with PKl, PKIR and PKDD [51] in comparison with those with 
TMl [50], NL3 [80], TW99 [61], DD-MEl [62], 



Interaction 


Psat.(fm~^) 


Eb [MeV] 


K [MeV] 


J [MeV] 


M*/M{n) 


M*/M(p) 


PKl 


0.1482 


-16.268 


282.644 


37.641 


0.6055 


0.6050 


PKIR 


0.1482 


-16.274 


283.674 


37.831 


0.6052 


0.6046 


PKDD 


0.1496 


-16.267 


262.181 


36.790 


0.5712 


0.5706 


NL3 


0.1483 


-16.250 


271.729 


37.416 


0.5950 


0.5950 


TMl 


0.1452 


-16.263 


281.161 


36.892 


0.6344 


0.6344 


TW99 


0.1530 


-16.247 


240.276 


32.767 


0.5549 


0.5549 


DD-MEl 


0.1520 


-16.201 


244.719 


33.065 


0.5780 


0.5780 



Although for the nonlinear self-coupling effective interactions, the density dependencies are only embodied 
in the Klein-Gordon equations, it is still worthwhile to obtain a quantitative understanding of the coupling 
constants. 

In Fig. ^ the density dependencies of the interaction strengthes for ga-{top), i^^^ (middle) and i^p (bottom) 
in nuclear matter as functions of the nuclear density are shown. The curves in the figures from top to bottom 
are labelled in the order of from left to right. The shadowed area corresponds to the empirical value of the 
saturation point in nuclear matter, i.e., Fermi momentum kp = 1.35 it 0.05 fm~^ or density p = 0.166 it 0.018 
fm~^. For the nonlinear effective interactions, the "equivalent" density dependence of the interaction strengthes 
for (7, a; and p is extracted from the meson field equations according to [51]: 

5(7 ~ 9(7 + [520-^ + 530-^] M, (56a) 
gu, ~ gcj -Csujq/Pv, (56b) 
9p ~ 9p - dsiplf/ps- (56c) 

For the cj-meson, TW-99 and DD-MEl exhibit quite different behaviors for compared with those of the other 
effective interactions in either magnitude or slope. In particular, the strengthes in TW-99 and DD-MEl for 
the density interval in Fig. |l]are almost two times larger than that of GL-97. Quite different results can also 
be seen at the empirical nuclear matter densities. For the w-meson, except for TW-99, DD-MEl, TMl, and 
GL-97, the strengthes in all the other effective interactions are density-independent. However, the strengthes 
are closer to each other at the empirical saturation density than those of the cr-meson although large differences 
can also be seen at low densities. For the p-meson which describes the isospin asymmetry, the strengthes for 
TW-99 and DD-MEl show strong density dependence in contrast with those of the other effective interactions; 
while those of PKl, PKIR and PKDD are just in between [51]. 

In Fig. |SJ the behavior of the binding energy per particle E/A as a function of the baryonic density p is 
shown for symmetric nuclear matter (left) and neutron matter (right). It is seen that all the density-dependent 
meson-nucleon coupling effective interactions give softer results than the nonlinear self-coupling ones , especially 
for neutron matter. The behaviors predicted by PKl, PKIR are much softer than that by NL3 and a little 
harder than that by TMl. The results from PKDD are slightly softer than that from DD-MEl and much 
harder than that from TW99 at high densities. All these behaviors can be explained in the density-dependent 
meson-nucleon coupling framework. 

As what have been mentioned in expressions H56|) , the meson-nucleon coupling constants in the nonlinear self- 
coupling of mesons can be expressed as some kind of density-dependence. Fig. 0] shows the density-dependence 
of the coupling constants for the nonlinear self-coupling effective interactions and for the density-dependent 
version, where almost all the density-dependent coupling constants decrease with increasing density except for 
of NL3, NLSH, NLl, which has strong self-couplings. On the other hand, the coupling constants and g^j 
of TMl, which has relatively weak a self-coupling {g^ = 0.6183) and strong uj self-coupling (cs = 71.3075), are 
smaller than the others, which means that TMl provides relatively weaker scalar and vector potentials. This is 
the reason why TMl presents the softer behavior than the other nonlinear self-coupling effective interactions. 
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Figure 4: Density dependence of the interaction strengths for <^o-(top), (middle) and i^p (bottom) in nuclear 
matter as functions of the nuclear density. The shadowed area corresponds to the empirical value of saturation 
point in nuclear matter (Fermi momentum kp = 1.35 ib 0.05 fm~^ or density p = 0.166 it 0.018 fm~^). Taken 
from Ref. [51]. 

In Fig. TW99 predicts the softest results because of its relatively small g^^ as compared with DD-MEl, 
PKDD and NL3, and large as compared with PKl, PKIR and TMl in Fig. ^ As it is known, the repulsive 
potential would be dominant at high densities. In Fig. [3 NL3 gives the hardest results because of its constant 
and large g^^ even though its <t-N coupling constant g^ increases with the density. For PKl, PKIR and 
PKDD, which present the mid soft behaviors in Fig El the coupling constants also lie between the largest and 
the smallest in Fig. |1J There exists a significant difference between symmetric nuclear matter and neutron 
matter. The density-dependent effective interactions PKDD and DD-MEl present similar trends as PKl, 
PKIR and TMl in symmetric nuclear matter but much softer in neutron matter, which may be interpreted 
by the density-dependence of gp. For the effective interaction PKIR, the density-dependence of the gp is fairly 
weak as compared with that of the density-dependent meson-nucleon coupling effective interactions. It can be 
explained by a very weak p-field, which generates neutron-proton asymmetry field. The behavior of gp with 
respect to the neutron-proton ratio is shown in FigEl As expected, the behavior is symmetric with respect to 
\n{N/Z) and the density-dependence becomes more remarkable with the increase of the baryonic density and 
the neutron-proton asymmetry. 

3 Relativistic continuum Hartree-Bogoliubov theory 

Pairing correlations are due to the short range part of the nucleon-nucleon interaction and play important 
roles in open shell nuclei. A simple and commonly used method of dealing with pairing correlations is the BCS 
approach under the constant gap approximation. The BCS method can be combined easily with the relativistic 
mean field theory as shown in Refs. [21,78]. However, the conventional BCS method is not justified for exotic 
nuclei because it could not include properly the contribution of continuum states. Bogoliubov transformation 
is the generalization of the BCS scheme. By quantizing the meson field and making Gorkov factorization, the 
relativistic Hartree Bogoliubov (RHB) formalism was derived [95]. In order to describe the exotic nuclei, the 
relativistic continuum Hartree Bogoliubov (RCHB) theory was developed, in which the continuum states are 
discretized and the RHB equations are solved in the coordinate space [36,42]. 

In this section, the pairing correlations, the conventional BCS approach, the continuum states, the resonant 
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Figure 5: The energy per particle EjA in symmetric nuclear matter (left) and neutron matter (right) as 
functions of the baryon density p. The shaded area indicates the empirical saturation region. The filled 
diamond presents the data taken from Ref. [94] as comparison. 
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Figure 6: The density-dependence of gp for PKIR with respect to the neutron-proton ratio N/Z. Taken from 
Ref. [51]. 
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BCS method and several methods for obtaining single particle resonant states are briefly reviewed. Finally, 
the detailed formalism of the RCHB theory together with the discussion on the effective pairing interactions 
are given. 



3.1 Pairing correlations 
3.1.1 BCS approximation 

For Fermion systems like atomic nuclei, Kramers degeneracy ensures the existence of pairs of degenerate, 
mutually time-reversal conjugate states, which could be coupled strongly by a short-range force, namely a 
pp-channel in the framework of mean field. Physically, such pairing correlations lead to the superfiuidity. In 
Hartree approximation, the ground state properties could be described by filling the single particle levels from 
the bottom up to the Fermi level. The occupation probability of each single particle level is either zero or 
one. Due to pairing interaction, pairs of nucleons are scattered from the levels below the Fermi level to those 
above. Thus for the open shell nuclei, one has to deal with the occupation probabilities ranging from zero to 
one. Mathematically, this could be treated by introducing the concept of quasi-particles. For stable nuclei, the 
pairing gap can be extracted from the experimental odd-even mass difference and the BCS approximation is 
very simple and useful. For M-single particle orbitals in the Hartree approximation, for each particle orbital k 
(fc = 1, • • • , M), one can introduce the quasi-particle energy: 

Ek=[{ek-Xf + A'f\ (57) 
and the corresponding occupation probability: 

efe - A 



2 1 



Eh 



(58) 

M 



where A is the Fermi level determined by the particle number N = f ^ and A the pairing gap determined 



from the experimental odd-even mass difference 



k=l 



A = i [E{N + 2) - E{N + 1) - {E{N + 1) - E(iV))] , (59) 

or deduced from certain pairing Hamiltonian //pair when the experimental binding energies are unknown. For 
monopole pairing interaction i/pair = G a'^a'^a^,ak', there lies the pairing gap equation with classical 



kk'>0 



variation: 



^ = G^UkVk = ^^^, (60) 

fc>0 fc>0 

where Uk satisfies = 1 — vf,. For fixed values of the gap parameters A or pairing interaction strengthes G 
for neutrons and protons, the self-consistent solution is obtained by iteration. 

The BCS method can be easily implemented into the relativistic mean field description with the single 
particle orbitals in the relativistic Hartree level. However, it is only valid for bound states. For the extremely 
neutron-rich (or proton-rich) nuclei near the drip lines, the continuum states begin to contribute as schematically 
shown in Fig. The conventional BCS approach will involve unphysical states and lead to the divergence. 
For these nuclei, one must either investigate the detailed properties of continuum states and include the 
coupling between the bound state and the continuum by extending the BCS method to resonant BCS theory 
(see Section I3.2.2|) or generalize the BCS model by unifying the Hartree (or Hartree-Fock) scheme and the 
Bogoliubov transformation (see Section [3. 



3.2 Continuum states 

3.2.1 Single particle resonant states 

For stable open shell nuclei, the pairing gaps are of the order of MeV [16]. There are some evidences that 
the pairing gaps will increase towards the drip lines. For nuclei close to drip lines, the one or two particle 
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Figure 7: The schematic picture for the difference in the pairing correlation of stable and drip line nuclei. 
Normally the last filled nucleon in stable nuclei has 8 MeV binding, while it is near the threshold in drip line 
nuclei. The pairing correlation in exotic nuclei provides the possibility to scatter valence nucleons back and 
forth in the continuum. 

separation energies are usually very small. For example, the one proton separation energies Sp are 0.138 MeV, 
0.60 MeV and 0.29 MeV in ^B, ^^N and ^^Cl; the two neutron separation energies S2n are 0.30 MeV and 1.33 
MeV in ^^Li and ^^Be [89]. A small particle separation energy means that the Fermi level is very close to the 
threshold. Therefore, pairing correlations can scatter valence nucleons between bound and continuum states as 
seen from Fig. [3 An unphysical continuum state distributes over the whole space while a physical resonance 
state with positive energy stays inside the nucleus for a long time and has better asymptotic behavior. For 
the cases where the pairing correlations involve the continuum, such as exotic nuclei, it has been demonstrated 
that the largest contribution to the pairing correlations is expected to come from the resonant continuum 
part [36,37,41,42,96,97]. Therefore it's very important to investigate single particle resonant states in exotic 
nuclei. 

Various techniques have been developed to study resonant states in the continuum such as the R-matrix 
theory [98] and the extended R-matrix theory [99], the K-matrix theory [100] and the conventional scattering 
theory [101]. By discretizing the continuum, the contribution of the resonant states can be self-consistently 
taken into account via a Bogoliubov transformation in coordinate space [36,37,42]. The bound-state-type 
methods have also been developed, including the real stabilization method [102], the complex scaling method 
(CSM) [103], and the analytic continuation in the coupling constant (ACCC) method [104]. Within the 
framework of the non-relativistic Schrodinger equation, the ACCC method has been applied to investigate 
the energies and widths for resonant states in light nuclei combined with few-body methods [105, 106], and 
to study single-particle resonant states in spherical and deformed nuclei by solving the Schrodinger equation 
with Woods-Saxon potentials [107]. Recently the ACCC method and phase shift scattering method have been 
combined with the RMF model. 

The philosophy of the ACCC method is that a resonant state becomes a bound one if one increases the 
attractive potential. Then the energy, width, and wave function of the resonant state can be obtained by an 
analytic continuation carried out via a Fade approximant (PA) from the bound-state solutions. Compared with 
other bound-state-type methods, the ACCC approach is very effective and numerically quite simple because 
many methods available for bound-state problems can be used and the PA for analytic continuation can be 
implemented in easily. 

By solving the spherical relativistic Hartree equations in a meshed box of size Rq self-consistently, one 
can calculate the ground state properties of a nucleus. The vector potential V{r) and the scalar potential 
S{r), energies and wave functions for bound states are also obtained. By increasing the attractive potential 
{V{r) + S{r)) to X{V{r) + S{r)) , a resonant state will be lowered and becomes bound if the coupling constant 
A is large enough. Near the branch point Aq, defined by the scattering threshold /c(Ao) = [104], the wave 
number ^(A) behaves as 




(61) 
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These properties suggest an analytic continuation of the wave number k in the complex A plane from the 
bound-state region into the resonance region by Pade approximant of the second kind (PAII) [104] 



k{x) « fci^'^i(x) = (62) 

where x = ^JX — and cq, ci, . . . , c^, di, ^2, . . . , are the coefficients of PA. These coefficients can be de- 
termined by a set of reference points and k{xi) obtained from the Dirac equation with Aj > Aq, « = 
1, 2, L -|- -|- 1. With the complex wave number A;(A = 1) = fc^ -|- ifcj, the resonance energy E and the width 

V can be extracted from the relation e = E — i— {E, F G M) and fc^ = — M^, i.e., 



(M2 + A;2 - A;2)2 + 4A:2A:2 + (^2 + A;2 - kf) 
E = \ ^ M, 

r = ^2^{M^ + k^-kfy + 4k^kf-2{M^ + k^-kf). (63) 

Based on the Schrodinger and Dirac equations, the stability and convergence of the energies and widths for 
single-particle resonant states with square-well, harmonic-oscillator and Woods-Saxon potentials have been 
investigated and their dependence on the coupling constant interval and the order of the PA have been dis- 
cussed [108, 109]. Combined with the RMF theory, the ACCC method was employed to study energies and 
widths for resonant states in stable nuclei ^^O and ^^Ca in Ref. [110] and to explore the single-particle resonant 
states, including resonance parameters and wave functions, in exotic nuclei in Ref. [111]. 

The continuation in the coupling constant can be replaced by the continuation in k plane along the A;(A) 
trajectory determined by Eq. 1)62(1 to the point kfi corresponding to the wave number for Gamow state, i.e., 
k^ = k^^'^^iX = 1) [104]. Similarly, the wave function ip{kji,r) for a resonant state can be obtained by an 
analytic continuation of the bound-state wave function ip{ki,r) in the complex k plane. One can also prove 
that the wave function ip{k, r) is an analytic function of the wave number k in the inner region r < Rq where 
the Jost function analyticity dominates [104]. Therefore, the technique, which has been adopted to find the 
complex resonance energy, has been used to determine the resonance wave function ip{kji,r). Firstly, the PA 
are constructed to define the resonance wave function at any point r in the inner region (r < Rq) [104] 

[L,N](, , PL{k,r) ^ ao(r)+ai(r)A; + 02(r)fc2 + ... + fl^(r)fc^ 

^ '""^ QN{k,r) l + b^(r)k + b2{r)k^ + ... + bN{r)k^ ' ^ ' 

where the coefficients ai{r) (i = 0, 1, . . . , L) and bj{r) (j = 1, 2, . . . , A^) are dependent on r. These coefficients 
can be determined by a set of reference points ki and ip{ki,r) obtained from the Dirac equation with Aj > 
Ao, {i = 1, 2, A^-|- 1). The resonance wave function </?(/c_R, r) = ip^^''^\kii,r) (r < Rq) can be extrapolated 
in this way. 

The Dirac equation can be rewritten as two decoupled Schrodinger-like equations for the upper and the 
lower components respectively [32]. For neutrons, the Schrodinger-like equation for the upper component reads 

- ^ii^Gl^ (r) + {Ef - M')G^(r) = 0, (65) 
in the outer region where V{r) ~ and S{r) ~ 0. Its solution is the well known Riccati-Hankel function 

h^{z) = h^{z) ± ijf,{z), z = kr, (66) 

where juiz) = zjt^{z) is the usual regular solution, i.e. Riccati-Bessel function, and ht^{z) = zn^iz) is the 
irregular solution, i.e. Riccati-Neumann function. The outer wave function is matched to the inner wave 
function a,t r = Vm < Rq 

^^^'^\kR,rm) = C{kR) Vj.ikRrm) + D{kR)n^{kRrm) \ , (67) 
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where C{kR) is the coefficient for matching and D{kji) = tand^^kji) with S^^kn) the phase shift. It has been 
found that S^ikji) is almost a constant when is large enough. Given the upper component, the lower 
component Fi^(r) can be calculated from the relationship between the upper and the lower components and the 
resonance wave function is finally normalized according to the Zel'dovich procedures [104]. The wave function 
for the neutron resonant state in ^'^Ca with different matching points is given in Fig. [H] where one finds 

the convergence of the wave function with respect to the matching point when changes from 14 fm to 
18 fm. 




5 10 15 20 25 30 35 40 45 50 
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Figure 8: Real parts of the upper and lower components of the radial wave function for the neutron 
state in ^'^Ca calculated by the ACCC method with different matching points r^- Taken from Ref. [111]. 

The resonant states can be also investigated in the scattering phase shift method [41, 97, 112] where the 
RMF equations are solved with the scattering- type boundary conditions. At large distances, where both the 
scalar and the vector potentials are zero, the RMF radial equations can be written in the form: 

(PG , 9 k;(k + 1),^ , , 

^ + - ^7^)G = 0' (68) 

" ' + (69) 



E + dr r 

where = -E^ — Af^. These equations are suited for fixing the scattering- type boundary conditions for the 
continuum spectrum. They are given by: 

G = Car[cos{5)ji{ar) — sm{5)ni{ar)], (70) 



Ga^r 

F = ^ ^ j^ [cosi6)ji^iiar) - sm{6)ni^i{ar)], (71) 

where j; and rii are the Bessel and Neumann functions and 6 is the phase shift associated to the relativistic 
mean field. The constant G is fixed by the normalisation condition of the scattering wave functions and the 
phase shift 5 is calculated from the matching conditions. In the vicinity of an isolated resonance the derivative 
of the phase shift has a Breit-Wigner form, i.e. 

dm ^ r/2 

dE {Er-Ey + ryA' ^ ' 

from which one estimates the energy and the width of the resonance. In the vicinity of a resonance the radial 
wave functions of the scattering states have a large localisation inside the nucleus. Close to a resonance the 
energy dependence of both components of the Dirac wave functions can be factorized approximatively by a 
unique energy dependent function [113]. As in the non-relativistic case [114], this energy dependent factor is 
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Figure 9: Energies and widths for the neutron states J^3p3/2, z^3pi/2, 2^2/7/2, 2^2/5/2, z^l/i9/2, and J^lii3/2 in 
^^^Zr. Sohd circles represent the results of the ACCC method, while open circles denote the results of the 
scattering method. Taken from Ref. [111]. 



the square root of the Breit-Wigner function written above, or, equivalently, the square root of the derivative 
of the phase shift. 

The energies and widths for the neutron resonant states J^3p3/2, t^'ipi/2j ^'^17/2^ ^2/5/2, i'lhg/21 and 1^1^13/2 
in "^^^Zr are shown in Fig. |^ The results of the ACCC and the scattering methods are in good agreement with 
each other for most of these states. From both methods, 2^2/5/2 and 1^1*13/2 have large widths, while 2^2/7/2 
and Z/I/19/2 are very narrow. The width for i'3p^/2 from the scattering method is slightly larger than that from 
the ACCC calculation. For the resonant state z^3pi/2) neither the energy nor the width can be extracted from 
the scattering calculation. From a quantum mechanical point of view, these single-particle resonant states are 
quasi-stationary ones captured by centrifugal barriers. The decay width for a resonant state can be roughly 
explained by the penetration through the barrier. For those states with the same I, i.e., the same centrifugal 
barrier, the higher state has a larger width, e.g., for the two / states 2^2/5/2 and 2^2/7/2- Although 1/I/19/2 is 
higher than 1/2/5/2, its width is smaller because of higher centrifugal barrier. A similar argument also holds for 
a much higher but narrow state z^lix3/2- 

3.2.2 BCS approximation with resonant states 

The conventional BCS method is incapable of describing weakly bound nuclei due to the oscillating asymptotic 
behavior of single particle wave functions in the continuum [37] . Recently, many authors demonstrate that by 
picking up only those low-lying resonant states in the continuum, the BCS method can be successfully applied 
to weakly bound nuclei [41,97,115-117]. This is because the resonant states are well localized inside the nucleus 
and there is a large region outside the nucleus where the resonant wave functions have values close to zero 
before they start oscillating [97]. Therefore, the wave functions of the resonant states can be taken to be zero 
beyond a cutoff radius [97]. Meanwhile a realistic pairing potential, such as a zero-range (5-force, prefer to pick 
up the resonant states, whose wave functions have a large overlap with those of bound states below the Fermi 
surface, rather than the continuum [115,117]. In such a way, the HF+rBCS and RMF+rBCS results do not 
depend sensitively on the cutoff radius [97]. 

The BCS method can be easily extended to study the width effects of resonant states. The extended BCS 
equations for a general (finite range) pairing interaction including the contribution of the resonant continuum, 
referred as the resonant-BCS (rBCS) equations [41,97], are : 

= E %J^J-^J- + E / gu{eW{eH{e)de , (73) 
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Figure 10: The rms neutron radii (left panel) and the two-neutron separation energies (right panel) of even Zr 
isotopes as a function of the mass number A. Taken from Ref. [41]. 



+ Yl Ke,i7F^ye^,177:; / 9u' {e)uu' {e'W {e)de' , 

I Jit 



N 



gl{e)vl{e)de . 



(74) 



(75) 



Here Aj are the gaps for the bound states and /S.^ are the averaged gaps for the resonant states. The quantity 

g'^ie) = ^ is the total level density and by is the phase shift of the state with the angular momentum 

vr ae 

(lyju)- The factor 5^(e) can take into account the effect of the width and it is approximately a delta function for 
a very narrow resonance. The interaction matrix elements are calculated with the scattering wave functions at 
resonance energies and normalized inside a volume where the pairing interaction is active (see last paragraph). 

The BCS equations in the conventional RMF+BCS formalism can be easily substituted by these rBCS 
equations (|7.S|74|75|) . To make this implementation, one can solve the RMF equations with the scattering- 
type boundary conditions or any other methods like the ACCC method described in the last subsection and 
find the energies, widths and wave functions for the resonant states. This approximation scheme, called the 
RMF+rBCS method, has been applied to study neutron-rich Zr isotopes [41]. It was demonstrated that the 
sudden increase of neutron radii close to the neutron drip line depends on a few resonant states close to the 
continuum threshold. Including these resonant states, the RMF-BCS calculations give practically the same 
neutron radii and neutron separation energies as the RCHB calculations (see Fig|lU|). 



3.3 Relativistic continuum Hartree-Bogoliubov theory 
3.3.1 Bogoliubov transformation and continuum spectra 

The contribution of the resonant states can be self-consistently taken into account in the Hartree-Fock Bo- 
goliubov (HFB) theory solved in the coordinate space [37]. The advantage of the general HFB theory is that 
the variation method based on quasi-particle transformation unifies the self-consistent description of nuclear 
orbitals and the BCS pairing theory into a single variation theory. In Fig. 1111 schematic pictures for different 
bases in the HF, HFB and the canonical basis have been given. For a given many-body system, one can choose 
any M-dimension orthogonal and complete basis to diagonalize the Hamiltonian. That is shown in the first 
column of Fig. 1111 If only HF approach are considered, the Hamiltonian is diagonal in the HF basis. We get a 
M-dimension basis, with a Fermi surface which divides the full open orbitals and full occupied orbitals. They 
are given in the second column. After the Bogoliubov transformation, we are working in the quasi-particle 
basis, which is enlarged into 2M-dimension. The quasi-particle orbitals are reflection symmetric with respect 
to the Fermi surface, as shown in the third column. They must be transformed into the physical space — 
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Figure 11: The schematic pictures of different basis. The first column is an arbitrary M-dimension orthogonal 
and complete basis to diagonalize the system Hamiltonian. The second column is the single particle levels in 
the HF space. The third column is quasi-particle levels in the quasi-particle HFB space which is enlarged into 
2M-dimension. The quasi-particle orbitals has a reflection symmetry with respect to the Fermi surface. The 
physical orbital in the canonical basis is given in the forth column. 




canonical basis [35,42], in order to have a similar explanation as the HF space in the second column. For the 
Bogoliubov transformation involving the continuum, as the Fermi energy is negative for a bound nucleus, only 
those resonant states in the canonical basis with proper asymptotic behavior could be picked up. In the forth 
column, a schematic picture for the orbitals in the canonical basis is given, now the orbitals could be partially 
occupied, which is another essential difference between the HF basis and the canonical basis. 

3.3.2 The formalism for the RCHB theory 

In the BCS approach, the gap equation takes into account only pairing between time-reversed continuum 
states. A more general pairing between continuum states at neighboring energies can be taken care of by a 
continuum Bogoliubov approach. The improvement to the BCS approximation is to introduce the concept 
of quasi-particle by the Bogoliubov transformation. Instead of introducing a mapping factor, the Bogoliubov 
transformation transforms the equation of motion into quasi-particle space. In the following we shall consider 
the transformation from the single particle basis to the quasi-particle basis. The single particle operators in 
coordinate space ar,ap, are connected to the operators in the quasi-particle basis. Pi, pj, i = 1,M, by the 
following transformation: 



where i and r represents the row and column indices respectively. The basis of the quasi-particle states, ■0c/,j(r) 
and V'yjilr)) are defined via the single particle coordinate basis of 5- function as: 

fa^'Sir-r') (^^A . (77) 



>y,i(r)y J V^i.r'/ 

In spherical case, the Dirac spinor wave functions ij)u^i{r) and ■0y,j(r) are similar to Ea. H2U() and satisfy the 
following normalization relation: 

' d^rW^W = 1. (78) 



The variation method based on quasi-particle transformation can unify the self-consistent description of nu- 
clear orbitals, as given by the HF approach, and the BCS pairing theory into a single variation theory. As the 
equations of motion are self-consistently solved in the quasi-particle space, the convergence is guaranteed auto- 
matically. The unphysical continuum are excluded and the contribution from the resonance states with positive 
energies can be taken into account. From now on, the concept of continuum in Bogoliubov transformation is 
just the resonance states with positive energy as the unphysical continuum are excluded already. 



26 



Following the standard procedure of Bogoliubov transformation, a Dirac Hartree-Bogoliubov equation could 
be derived and the unified description of the mean field and pairing correlation in nuclei could be achieved. 
Using Green's function techniques it has been shown in Refs. [21,95] how one can derive a relativistic Hartree- 
Fock-Bogoliubov theory: after a full quantization of the system the mesonic degrees of freedom are eliminated 
and, in full analogy to the non-relativistic case, the higher order Green's functions are factorized in the sense 
of Gorkov [118]. Finally, neglecting retardation effects and the Fock term, as it is mostly done in relativistic 
mean field theory, one ends up with RHB equations as the following: 

^ ](*") =e(*"), (79) 

\^ -A* -h* + x J Ki^vJ Ki'vj 

where 

/i(r, r') = [a • p + V{v) + /3(M + 5(r))] (80) 
is the Dirac Hamiltonian and the Fock term has been neglected. The pairing potential is 

Afcfc,(r,r') = -y"(i\iy"a!M^V;,,^^,(rr';rir;)^^^,(ri,r;), (81) 

kk' 

where k, k' , k, or k' represent the other quantum numbers which together with the coordinate specify the single 
particle states. It is obtained from the pairing tensor, 

/^fcfc'(r,r') = {\ak,iak',i'\) = ^ ^/;^'*(r)*Vy •*(r'), (82) 

E^>0 

and one-meson exchange interaction V^^, (rr'; rir'^^) in the pp-channel. The nuclear density is given as 

p(r,r')= Yl i^'v\rn'/ir'). (83) 

k,Ei>0 

In the RHB theory, the ground state |^) of the even particle system is defined as the vacuum with respect 
to the quasi-particle: = for all i, i.e., |^) = rij/^*!")) where |— ) is the bare vacuum. For an odd 

particle system, the ground state can be correspondingly written as: = /?j Yiij^j where j is the level 

which is blocked. The exchange of the quasiparticle creation operator /?j with the corresponding annihilation 
operator Pj means the replacement of the column j in the U and V matrices by the corresponding column in 
the matrices V*, U* [35]. 

As the one-meson exchange interactions are not able to reproduce even in a semi-quantitative way the 
proper pairing in the realistic nuclear many-body problem, the interaction V in Eq. ()81() is replaced by the 
phenomenological interaction which has been generally used in the pairing channel of the conventional HFB 
theory. As in Ref. [36] V used for the pairing potential ()81() is either the density-dependent two-body force of 
zero range: 

V{ri,r2) = VoS{ri - ra)^ [1 - cricra] (l - ^) (84) 
with the interaction strength Vq and the nuclear matter density po or Gogny-type finite range force: 

V{ri,r2) = Ye^^^'-^^^/i^'^^iWi + B^P" -HiP^ -AUP^P^) (85) 

i=l,2 

with the parameter /ij, Wi, Bi, Hi and Mj (i = 1,2) as the finite range part of the Gogny force DIS [119]. 

The RHB equations can be solved in different bases as well. In many applications an expansion of the wave 
functions in an appropriate harmonic-oscillator basis of spherical or axial symmetry provides a satisfactory 
level of accuracy. For example, in Ref. [120] the RHB equations were solved by expanding the nucleon spinors 
Uk{r) and 14 (r) and the meson fields in a basis of spherical harmonic oscillators. However, for nuclei near the 
drip lines the expansion in the localized oscillator basis presents only a poor approximation for the continuum 
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states. Oscillator expansions produce densities which decrease too steeply in the asymptotic region at large 
distance from the center of the nucleus. The calculated rms radii cannot reproduce the experimental values, 
especially for halo nuclei. The transformed harmonic oscillator (THO) basis has also been employed in the 
solution of the RHB equations in configurational space [64]. Solving the RHB equations in Woods-Saxon basis 
is another alternative. 

In order to describe the coupling between bound and continuum states more exactly, the RHB equations 
and the equations for the meson fields should be solved in coordinate space [36,42, 121] and one arrives at the 
RCHB theory. When spherical symmetry is imposed, the wave function can be conveniently written as 



u 



YUO,(t>)Xt{t), ij'v 



\ 



YUO,<P)xt{t)- 



(86) 



The above equation (|79jl only depends on radial coordinates and can be derived as the following integro- 
differential equations: 



r dGujr) 

dr 
dFujr) 

dr 
dGvjr) 

dr 
dFyjr) 

dr 



+ '^Gu{r)-{E + \-V{r) + S{r))Fu{r)+r j r'dr' A{r,r')Fvir') 

- '^Fu{r) + {E + X- V{r) - S{r))Gu{r) +r j r'dr'A{r, r')Gy (r') 
+ '^Gv{r) + {E-X + V{r) - S{r))Fvir) +r j r'dr'A{r, 

- '^Fvir) -{E-X + V{r) + 5(r))Gy (r) +r j r'dr'A{r, r')Guir') 



0, 
0, 
0, 
0, 



(87) 



where the nucleon mass has been included in the scalar potential S{r). Instead of solving Eqs. (|2H) and (|2r 
self-consistently for the RMF case, now one has to solve Eqs. ()87() and H23() self-consistently for the RCHB 
case. The densities are calculated as 



47rr^Ps(r) 
47rr^/9^(r) 
47rr^/93(r) 
Airr'^Pcir) 



E.(|GV(r)p-|F(.(r)|2), 

Eii\Gv{rr + mr)n 

E.r3(|GV(r)|2 + |F^(r)|2), 
E^hi^-n){\GUr)\' + \Fl.{r)\^) 



where the summations are over all quasi particle states. From the densities given above, one can calculate the 
rms radii and charge radius of the nucleus. The total binding energy is calculated as 
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E. 



nucleon 



+ E^ + E^ + E„ + E, + E, 



CM, 



where except for -Enucieom the other terms are the same as those in the spherical RH theory. The energy of 
nucleons -Enucieon is calculated as 



(90) 



and E, 



pair 



— |TrAK. For the details of the transformation from the quasi particle space to the canonical 
basis, the reader is referred to Ref. [42]. 

To solve RCHB equations (|87|) . one has to discretize the RCHB continuum. The discretization can be done 
by solving the RCHB equation in a spherical box of radius R, i.e., by imposing the boundary conditions 



i^uiE, |r| =R,s)= i^viE, |r| = R, s) 



0. 



(91) 



By increasing R, one can make the RCHB spectrum dense and better approximate the continuum, while 
increasing Emax allows one to take into account coupling to highly excited quasi-particle state. Both R and 
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Emax must be reasonably taken in such a way that the final results do not depend on the chosen R and Emax- 
Only those solutions with E > Q are necessary due to the special symmetry of RCHB equation [42]. If E is the 
solution of Eq. ((7^ . —E must be another solution, but they are not independent. Both solutions are connected 
as following: 

i^\j{-E,v)=^l,\,{E,v), i;\,{-E,r) = -i;ij{E,r). (92) 

So one can concentrate only on the solution with either positive energies or negative energies. 

Equations (jSTj), in the case of S-force Ea. (|5H) . is reduced to ordinary coupled differential equations and 
can be solved with shooting method combined with Runge-Kutta algorithms, which is so far the most elegant 
method for coupled differential equations in coordinate representation. For each fixed value of the energy E 
the RCHB equations (|87() are solved in the following steps: 

1. Discretizing the whole space [0,cxd]; 

2. Choosing a value for the energy E in Eq. (|87|) with E > for the given spin-parity channel; 

3. Finding the proper boundary condition for ip{0), where represents Gv,Fv,Gu,Fu and integrated 
Eq. H87() outward from r = to proper matching point r = Rmatch by Runge-Kutta algorithms; 

4. Finding the proper boundary condition for ^/>(oo) and integrated Eq. (|87|) from r = oo inward to proper 
matching point r = Rmatch by Runge-Kutta algorithms; 

5. Requiring the 'ip{r) derived from the former two steps to be the same at r = Rmatch will lead to the 
correct energy E; 

6. Repeating the process till all the energies in the given spin-parity channel lying in [0, Ecut] have been 
found. 

Equations (|87j) for Gogny forces (|85|l are a set of four coupled integro-differential equations. They are discretized 
in the space and solved by the finite element method as the following: 

1. Discretizing the whole space [0, oo] with + 1 points at ri = 0, r2 = p, r-jv+i = N x p, where p = R/N. 

2. Assuming ^(r) at ri = 0, r2 = /?, ...,rAr+i = N x p are known, if p = R/N is small enough, ^/>(r) in the 
element (i — l)xp<r<ixp could be written as a simple linear combination of ip{ri) and ip{ri^i); 

3. Requiring Eqs. ()87() is valid in the element (i — 1) x p < r < i x p, which could be easily integrated 
analytically and lead to a simple relation for E, ip{ri) and ip{ri^i); 

4. Repeating the process for all the element, one then is lead to a 4 x (A^ + 1) linear algebra equations 
for E, Gv,Fv,Gu,Fij for r = ri,r2, ...,rAr+i, which is nothing but a generalized eigenvalue problem; 

5. Diagonalizing this generalized eigenvalue problem one can get all the energies and their wave function 
discretized at r = ri, r2, ttv+i- 

Normally, equations 1)8 7|) and H23() are solved in a self-consistent way by the shooting method or finite element 
methods with a step size of 0.1 fm using proper boundary conditions in a spherical box of radius R < 20 fm. 
The Fermi-surface A is determined self-consistently by the particle number [42]. The results do not depend on 
the box size for i? > 15 fm [43]. For (5-force (|84|) . the number of continuum levels is limited by a cut-off energy, 
which must be larger than the depth of the potentials. 



3.3.3 Pairing correlation: finite range vs. zero range 

In the vicinity of the drip line, as the contribution from the continuum is switched on, it is very challenging to 
have the proper single particle levels and treat the coupling of the continuum without any ambiguity. Therefore, 
a proper interaction for the coupling between the bound state and the continuum is also essential. In principle, 
the effective nucleon-nucleon interaction should be obtained by means of Bruckner renormalization which gives 
the correct interaction after modifying the free interaction for the effect of the nuclear medium. In practice. 
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Figure 12: Left: Two-neutron separation energies S2n of even Ni isotopes as a function of A^, including the 
experimental data (solid points), RCHB with (5-force (open circles), RCHB with Gogny force (triangles) ,HFB 
with SkP interaction (stars) and HFB with SIII interaction (pluses); Right: The neutron pairing energies for 
RCHB with Gogny force (triangles) and 5-force (open circles). The open circles connected by dashed line or 
solid line are respectively the pairing energies for RCHB with (5-force by fitting that of RCHB with Gogny force 
at ^^Ni or ^^Ni. Taken from [122]. 



however, the effective interactions are approximated by some phenomenological forces, such as Skyrme-type 
(5-force or finite range Gogny force. Delta force is relatively simple for numerical calculation and has a realistic 
density-dependent behavior, but it allows a coupling to the very highly excited states. Normally an energy 
cutoff has to be introduced and the interaction strength has to be properly justified. The Gogny force has 
a better treatment for the coupling to the highly excited states, but it involves more sophisticated numerical 
techniques. 

The Skyrme- and Gogny-type pairing forces have been discussed before within the conventional HFB for- 
malism, but the HFB equations for Gogny force are solved on a harmonic basis, which is valid only for stable 
nuclei. In Ref. [122], the proper form of the pairing interaction near the drip line is discussed in the framework 
of the RCHB theory. The chain of even-even nickel isotopes ranging from the proton drip line to the neutron 
drip line are taken as examples. The pairing correlations are taken into account by either a density-dependent 
force of zero range or the finite range Gogny force. Through the comparisons of the two neutron separation 
energies S'2n, the neutron, proton, and matter rms radii, excellent agreements have been found between the 
calculations with both interactions and the empirical values. In Fig. 1121 for example, the two neutron separa- 
tion energies S2n of even Ni isotopes are shown as a function of the neutron number N from the proton drip 
line to the neutron drip line, including the experimental data (solid points), RCHB with (5-force (open circles), 
RCHB with Gogny force (triangles), HFB with SkP interaction [38] (stars) and SIII interaction [123] (pluses). 
The RCHB results with (5-force and Gogny force are almost identical. They show a strong kink at N=28 and 
a less one at A^=50. The drip line nucleus is predicated at ^'^''Ni in both calculations. The empirical data is 
known only up to A^=50. Comparing with the available empirical data, the general trend and gradual decline 
of S2n have been well reproduced. 

The neutron pairing energies — ^TrA/t for the calculation of RCHB with Gogny force are given in Fig. [T^ 
as triangles. The pairing energies demonstrate the shell structure revealed by the S2n in Fig. 1121 They vanish 
at the closure shell and have a maximum values in the middle of two neighboring closure shells. A tendency of 
the enhancement of pairing energies appears for the neutron-rich nuclei. The pairing energies for RCHB with 
(5-force are given by open circles. The circles connected by the solid line and dashed line are calculations with 
Vq = 742 MeVxfm"^ (by fitting the pairing energy of the Gogny force at ^'^Ni) and 715 MeVxfm"'^ (by fitting 
at ^^Ni), respectively. The variation behavior of the pairing energies with is quite similar, although Vq fitted 
at ^''Ni gives 2 MeV more for stable nuclei and Vq fitted at ^^Ni gives 4 MeV less for neutron rich nuclei. But 
these differences have little influence on the two-neutron separation energies and the rms radii, as one observes 
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The occupation Probability of Ni 




Figure 13: The occupation probabilities in the canonical basis for ^"Ni as a function of the single particle energy 
around the threshold, for RCHB with 5-force (dashed line) and Gogny force (solid points). Taken from [122]. 



in Fig. I12( for example. This clearly demonstrated the reliability of the result against the pairing strength. 

Another interesting feature of the exotic nuclei is the contribution from the continuum. Definitely it is 
very interesting to see how different the contribution from the continuum is for the calculations with Gogny 
force and 6-foicce. In Fig. the occupation probabilities of ^''Ni in the canonical basis is shown [35] for the 
neutron levels near the Fermi surface, i.e. in the interval — 10 < E^ p^ < 10 MeV. The occupation probabilities 
for RCHB (5-force are represented by the dashed line and RCHB with Gogny force by the solid points. For 
Vo = 715 MeVxfm"'^ or Vq = 742 MeVxfm"'^, the occupation probabilities are the same. Even the difference 
from the Gogny and the (5-force is too small to have any influence on the physical observables like S2n and rms 
radii. 

It has been shown that after proper renormalization (e.g., fixing the corresponding pairing energies), physical 
observables such as S2n and rms radii remain the same, even within a reasonable region of the interaction 
strength [122]. 



4 Exotic nuclei 



The study of exotic nuclei has attracted world wide attention due to their large N/Z ratios (isospin) and 
interesting properties such as halo and skin. It is also important for other fields including astrophysics, e.g., 
the properties of exotic nuclei are essential to understand the nucleosynthesis in the r process. Since the first 
case of halo in an exotic nucleus ^^Li was observed with RIB in 1985 [13], more and more exotic nuclei have been 
investigated with various modern experimental methods to understand this attractive phenomenon better [6-8, 
11,12]. For nuclei far from the /3-stability valley and with small nucleon separation energy, the valence nucleons 
in exotic nuclei extend over quite a wide space to form low density nuclear matter. Furthermore, the Fermi 
surface for exotic nuclei is usually very close to the continuum threshold. The valence nucleons could be easily 
scattered to the continuum states due to the pairing correlation. Thus, theories which can properly handle the 
pairing and continuum states are needed to describe the properties of exotic nuclei. One of such theories is the 
RCHB theory [36,42] introduced in the previous section. In this section, the RCHB theory will be applied to 
study the properties of exotic nuclei, e.g., the binding energies, particle separation energies, the radii and cross 
sections, the single particle levels, shell structure, the restoration of the pseudo-spin symmetry, the halo and 
giant halo, and halos in hyper nuclei, etc. 



4.1 Binding and separation energies 

The self-consistent microscopic RCHB theory has been used to study the ground state properties systematically 
such as binding energy, separation energy, radius, density distribution, single particle spectrum and shell 
structure etc. for a large number of isotope chains from the proton drip line to the neutron drip line throughout 
the periodic table including Li, C, N, O, F, Na, Ca, Ni, Zr, Sn and Pb [32,36,42-47,122,124-126]. In this 
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Figure 14: The two neutron separation energies S2n for the proton magic isotopes pfotted as a function of 
the neutron number in RCHB with NLSH. The experimental values are denoted by the solid symbols and the 
calculated ones by the open symbols. Taken from Ref. [126]. 



subsection, we will briefly present the description of the binding energies and separation energies for some 
typical spherical nuclei in RCHB theory. 

Ground state properties of all the even-even O, Ca, Ni, Zr, Sn, Pb isotopes ranging from the proton drip line 
to the neutron drip line were investigated with the RCHB theory [126]. The binding energy Ef^ calculated from 
the RCHB theory for these isotope chains and the corresponding data available [89] are compared. Generally 
speaking, the calculations reproduce the available data of binding energy quite well. The errors are usually 
within 3-4 MeV which is less than 1% of the experimental values. Since the deformation effects are not included 
in the RCHB calculations yet, for some nuclei with large deformation, larger discrepancies {AEf, ~ 10 MeV) 
of binding energy from experimental values are found, e.g., in some Zr isotopes. 

In Table the binding energies for Pb isotopes calculated with different effective interactions PKl, PKIR, 
PKDD, TMl, NL3, TW99 and DD-MEl are given in comparison with the experimental values [51]. Many 
effective interactions, particularly PKl, PKIR and PKDD provide good descriptions of the available data. All 
the effective interactions overestimate the binding energy in the beginning of the isotopic chain. However, from 
i9°Pb to 2i0pb, the effective interactions PKl, PKIR and PKDD give better descriptions than all the others. 
There exist however systematic deviations out of the neutron magic numbers, e.g., in ^^^Pb. For the effective 
interactions TMl, NL3, TW99 and DD-MEl, the deviations in ^i^pb are smaller but fairly large in ^ospb. 

The one or two nucleon separation energies or S2p{2n) ^'^^ quite sensitive quantities to test a microscopic 
theory. In Fig.^Jboth the theoretical and the available experimental S2n ai^e presented as a function of neutron 
number N for the O, Ca, Ni, Zr, Sn, Pb isotope chains, respectively. The good coincidence between experiment 
and calculation is clearly seen. Seen along the S2n versus N curve in Fig. I14| strong kinks can be clearly seen 
at iV = 8, iV = 20,28,40, N = 28,50, N = 50,82 and N = 126 for O, Ca, Ni, Zr, Sn, Pb isotope chains, 
respectively. All these numbers correspond to the neutron magic numbers, which come from the large gaps of 
single particle energy levels. 

For the proton magic even mass nuclei studied in Ref. [126], the neutron drip line nuclei are predicted as 
2*0, ■^^Ca, '^^Ni, ^3^Zr, ^^^Sn, and ^'^'^Pb, respectively. In Refs. [127,128], the last bound nucleus in neutron 
rich Ca isotopes is predicted as ™Ca with the HFB and Skyrme HF methods, respectively. The reason for that 
^^Ca is bound is caused by the halo effect, which also results in the disappearance of the normal = 50 magic 
number in the RCHB calculation as shown in Fig. 1141 where no kink appears at A^ = 50 for Ca isotopes. In 
stable nuclei, the A^ = 50 magic number is given by the big energy gap between the lg'9/2 and 3s, 2d levels. 
However in the vicinity of the drip line region of Ca isotopes, the single particle energy of 3,Si/2 would decrease 
with the halo effect and its zero centrifugal potential barrier, the large gap between 159/2 and 3si/2 disappears. 
Therefore, the nucleus ™Ca is no more a doubly magic nucleus and ''^Ca is bound. This disappearance of 
the A^ = 50 magic number at the neutron drip line is due to the halo property. This is different from the 
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Table 5: Total binding energies of Pb isotopes (in MeV) calculated with the nonlinear self-coupling effective 
interactions PKl, PKIR and the density-dependent meson-nucleon coupling effective interaction PKDD [51], 
in comparison with the experimental values [89] and the results of TMl, NL3, TW99 and DD-MEl. 
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Figure 15: The one neutron separation energies 5n for the nuclei ^^C/^ ^^N/^ ^^O and ^^F by RCHB 
theory(open circles) and their experimental counterparts (sohd circles). Taken from Ref. [46]. 



disappearance of the = 20 magic number which is due to deformation. Recently a new magic number 
= 16 has been discovered in neutron drip line light-nuclei region which is also considered as owing to the 
halo effect [14]. It can be expected that near drip line, more traditional magic numbers would disappear and 
new magic numbers would be found with the same mechanism. 

Systematic RCHB calculations with NLSH have been carried out for the C, N, O and F isotopes [46]. The 
one neutron separation energies predicted by RCHB and the data [89] for the nuclei ^^~^^C, ^'^"■^^N, 15-26q 
and i'^~25p g^j^g shown in Fig. as open and solid circles respectively. For carbon isotopes, the theoretical 
one neutron separation energies for ii-18.20,22q ^^^^ -^^ agreement with the data. The calculated Sn is less 
than (—0.003 MeV) for the odd-^d nucleus ^^C which is bound in experiment. While for the experimentally 
unbound nucleus ^^C, the predicted value of is positive. From the neutron-deficient side to the neutron 
drip line, excellent agreement has been achieved for the nitrogen isotopes. As almost all relativistic mean field 
approaches do, the RCHB calculations overestimate binding for 25,26q ^j^if.]^ g^j-g unstable in experiment. For 
fluorine isotopes, the overestimated in contrast with the underestimated one in °F. The 

last bound neutron-rich nucleus is predicted as ^''F. In general, the RCHB theory reproduces the data well 
considering that this is a microscopic and almost parameter free model. As discussed above, there are some 
discrepancies between the calculations and the empirical values for some of the studied isotopes. This may be 
due to deformation effects, which has been neglected in Ref. [46]. 



4.2 Radius, nuclear density distribution, and cross section 
4.2.1 Nuclear root mean square radii 

The nuclear radii are important basic physical quantities to describe atomic nuclei. The calculated neutron, 
proton, matter, and charge radii for all even Ca isotopes against mass number A are plotted in Fig. 1161 For 
nuclear charge radii Tc, the well-known parabolic behavior along ^'^"^'^Ca is not reproduced in RCHB calculation. 
It is mainly due to the improper spherical supposition to the real deformed ^■^•'^^'^^Ca nuclei. It can be clearly 
seen that the radii r^, r„, and Tc all increase with the mass number A. The mass radii as well as r„ 
increase much faster than rp and rc- In addition to the normal increase of r„ and with A, a slightly up-bend 
tendency occurs in the neutron drip line region from ^ = 62 up to A = 72. 

For light nuclei, the RCHB calculations predict a slight increase of proton or charge radii when the proton 
drip line is approached [46]. The neutron and proton rms radii predicted by RCHB for the nuclei ^''"^■^C, 
12-24]\^^14-26q g^^^ i6-25p ^^.^ given in Fig. 1171 The neutron radii for nuclei in each isotope chain increase 
steadily. While the corresponding proton radii remains almost constant for nuclei in each isotope chain except 
for proton rich ones. The increase of proton radii with neutron number decreasing in the proton rich side in 
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Figure 16: The rms radii for the proton, neutron, charge and matter distributions are pfotted as a function 
of the neutron number for the Ca isotope. The available experimental values are denoted by solid symbols. 
Taken from Ref. [126]. 
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Figure 17: The neutron and proton rms radii predicted by RCHB for the nuclei ^^C, ^^N,^^ ^^O and 
i6-25p^ Taken from Ref. [46]. 



each isotope chain is responsible for the slightly increase in the charge changing cross section of these nuclei on 
the target -^^C as discussed in Section 14.2.41 

In Fig. ^Jthe neutron radii r„ from the RCHB calculation for the even-even nuclei of the O, Ca, Ni, Zr, Sn 
and Pb isotope chains are ploted. The predicted rn curve using the simple empirical equation r„ = ro • N'^^^ 
with ro = 1.139 normalizing to ^'^^Pb is also represented in the figure. The simple formula for r„ agrees with 
the calculated neutron radii except for two anomalies. One anomaly appears in the Ca chain above N = 40, 
and the other in Zr chain above = 82. The increase of exotic Ni and Sn nuclei is not as rapid as that in Ca 
and Zr chains. The regions of the abnormal increases of the neutron radii are just the same as those for S2n- 
Both behaviors are connected with the formation of giant halo (see Section [4.41 for details). 

It is interesting to investigate the neutron number dependence of the proton radii (or charge radii) and 
the neutron radii for these proton-magic isotope chains. In Fig. 1191 the rms charge radii Tc obtained from the 
RCHB theory (open symbols) and the data available (solid symbols) for the even-even Ca, Ni, Zr, Sn and Pb 
isotopes are shown. The RCHB calculations reproduce the data very well (within 1.5%). For a given isotopic 
chain, an approximate linear dependence of the calculated rms charge radii Tc is clearly seen. However, the 
variation of Tc for a given isotopic chain deviates from the general accepted simple A^^^ law (denoted by dashed 
lines in Fig. 119(1 . which shows a strong isospin dependence of nuclear charge radii for nuclei with extreme N/Z 
ratio. Recently, these charge radii as well as the large amount of experimental data are investigated. A formula 
with the Z^/^ dependence and isospin dependence is proposed which describes well the charge radii all over the 
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nuclear isotope chart, see subsection 14.2.31 [129]. 
4.2.2 Neutron skin thickness 

The measurement of nucleon density distributions provides basic and important information in nuclear struc- 
ture. The proton distribution of stable nuclei can be determined accurately by elastic electron or muon scatter- 
ing. However, there is not comparable measurement of the neutron density distribution up to now. Although it 
is difficult to obtain the neutron density distribution, the difference in radii of the neutron and proton density 
distributions can be more easily determined [130]. Theoretically, the neutron skin thickness, r^ — rp could place 
important constraints on effective interactions used in nuclear models [131-133]. 

Many investigations have been done aiming at the description of the neutron skin thickness in the framework 
of the RMF model. With the pairing correlation treated by a constant BCS approach, the neutron rms 
radii of some closed and open shell nuclei are investigated with the spherical RMF model [134]. It has been 
demonstrated that with an increase in isospin of nuclei, the neutron rms radii in the relativistic mean-field 
approach with the parameter sets NLl and NL2 are larger than both the empirical data and the predictions of 
the Skyrme mean field [135]. This overestimation of the neutron skin thickness remains true for the effective 
interaction NL3, although the relativistic Hartree-Bogoliubov model with NL3-I-D1S give a good result of the 
neutron-skin thickness along the Sn isotopic chain compared to the available data [136,137]. 

The overestimation of the neutron skin thickness is attributed to a poor description of the isovector channel 
of the nonlinear effective interaction in Ref. [62] where the RHB model is extended to include density-dependent 
meson-nucleon couplings. A density-dependent effective interaction DD-MEl is proposed in which the improved 
isovector properties result in a better description of charge radii, and especially the calculated values of neutron 
thickness are in much better agreement with experimental data. As what is shown in Fig. 1191 the difference 
between the values calculated with NL3 and DD-MEl increases with the number of neutrons to about 0.1 
fm at = 82, but then it remains practically constant for A^ > 82. While both interactions reproduce the 
isotopic trend of the experimental data, NL3 obviously overestimates the neutron skin. The values calculated 
with DD-MEl, on the other hand, are in excellent agreement with the experimental data. This result presents 
a strong indication that the isovector channel of the effective interaction DD-MEl is correctly parametrized. 

In Ref. [133], the sensitivity of the neutron skin thickness S" = rn — rp in ^osp^j ^j^g addition of nucleon-o"-/? 
coupling corrections to a selection (PKl, NL3, S271, and Z271) of interactions in RMF model was investigated. 
With the higher order nucleon-o"-/) coupling, the PKl and NL3 effective interactions lead to a minimum value 
of the neutron skin thickness S = 0.16 fm and the S271 and Z271 effective interactions yield even smaller values 
oi S = 0.11 fm. The addition of these higher order correction terms to the conventional RMF models can be 
associated with some kind of multimeson exchange processes occurring in the inner (higher density) region of 
nuclei. In the inner region of a heavy nucleus such as ■^''^Pb, one could imagine that multimeson exchange 
processes introduced through nucleon-a-p coupling corrections could play an important role in modifying the 
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Figure 19: DD-MEl and NL3 predictions for the neutron skin thickness of Sn isotopes, compared with experi- 
mental data. Taken from Ref. [62]. 



short-range behavior of the nucleon-nucleon interaction, by softening the density dependence of the symmetry 
energy at high density (the inner part of nuclei) and thus leading to smaller values of neutron skin thickness S 
compared to previous effective interactions. 

4.2.3 Nuclear charge radii: isospin dependence and the Z^^'^ law 

Among all the size quantities describing nucleus, nuclear charge radius has been measured with very high 
precision by various techniques and methods experimentally. With its accuracy, the study of the nuclear 
charge radii is very important to understand not only the proton distribution inside the nucleus but also the 
exotic phenomena such as halo and skin observed in exotic nuclei. Particularly if one can get a simple and 
reliable formula for nuclear charge radii, it will be very useful to extract the decoupling of proton and neutron 
in the exotic nuclei and provide information for the effective nucleon-nucleon interaction widely used in all 
the nuclear models. The available charge radii data for ^ > 40 were examined and its global behavior was 
studied in Ref. [129] . Instead of the widely accepted A^^^ law, a new Z^^^ formula with isospin effect has been 
proposed. 

Based on the consideration of the nuclear saturation property, nuclear charge radius Rc is usually described 
by the A^^^ law [16,35] 

i?e = vaA^^^ (93) 

where A is the mass number and Rc = rc, with rc the rms charge radius. For very light nuclei, because 

of their small A and large fluctuation in charge distribution due to the shell effect with short period, it seems 
that the charge distribution radius as a bulk property has little meaning. A detailed analysis of charge radii 
data for j4 > 40 shows that is by no means a constant, but systematically decreases with A; i.e., « 1.31 
fm for light nuclei {A ~ 40) and ~ 1.20 fm for very heavy nuclei (see upper left panel in Fig. I20|) . This 
fact implies that some physics is missing in Eq. I|93p . 

Definite evidences of the violation of A^^^ law have been found in the measurements of isotope shift in rms 
charge radii [138-140]. In particular, 5{r'^)A+2,A values (associated with an addition of two neutrons) are often 
found to be considerably smaller compared to what is expected from the A^^^ law {6{r'^)A+2,A = ^ 
A typical example is that the observed charge radii of the calcium isotopes ^'^ ^''Ca remain almost the same 
(except a very little change induced by deformation or shell effect), though the mass number A has changed 
significantly. In contrast, there is also evidence that the observed 6{r'^)A+2,A values (associated with the 
addition of two protons) are often greater than what is expected from the A^^^ law (e.g., 5{r'^) for ^^Ti— ^"^Ca, 
50Ti_48ca, etc.). 

Along the /5-stability line, the ratio Z/A gradually decreases with A, i.e., for light nuclei Z/A ~ 1/2, and 
for the heaviest /3-stable nucleus ^SfU, {Z/A)^^^ « 0.7285. It turns out that (1/2)1/3/0-7285 « 1.09, which is 
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Figure 20: Comparison of descriptions of some formulas for nuclear charge radius Re- The parameters in the 
first four formulas are obtained through least square fitting to data on charge radii of nuclei with A > 40. 
Taken from Ref. [129]. 

very close to the ratio 1.30/1.20 shown in the upper left panel of Fig. 1201 A naive point of view is that the 
charge radius of a nucleus may be more directly related to its charge number Z, rather than its mass number 
A. Therefore, compared to the A^^^ law, a Z^^^ dependence for nuclear charge radii may be more reasonable, 

R, = TzZ^I^, (94) 

as noted in Refs. [138,139]. An analysis of the very limited data of charge radii then available showed that rz 
remains almost a constant, i.e., rz = 1.65(2) fm for A > 40. Similar Z^/^ dependent formula for nuclear charge 
radii was also proposed in Refs. [141, 142]. With the physics for such a simple Z^^^ law open, one notes here 
that the Z^^'^ law of Rc includes an isospin dependence compared to the A^^^ law as following: 

. ,..V,..,(|_^)-.^..3(,_^)- 



(95) 



One can get = rz/2^^'^ ~ 1.31 fm and 0.33 for the coefficient of the isospin term. Equation (|95j) is 
approximately equivalent to that obtained in Ref. [140], Rc = 1.25^^/^ (l — 0.2 ^^^ ) , by fitting the data. 
Obviously Eq. (jUlJ has merits that only one parameter is needed. 

The Z^^'^ dependence of nuclear charge radii was also used to modify the Coulomb energy term in the 
semi-empirical nuclear mass formula [143], and it was found that the agreement between the calculated and 
experimental results was improved. Moreover, the A~^^^ law for the nuclear giant (monopole, dipole and 
quadrupole) resonance energy (oc l/R) also could be improved, if the A~^^^ dependence is replaced by a Z~^/^ 
dependence [144]. 

In the past two decades, a vast amount of new experimental information on the electromagnetic structure 
of nuclear ground states of many nuclei has become available with high accuracy. The vast amount of improved 
experimental results have been investigated by using the A^^^ and Z^/^ dependence [129]. The measured 
charge radii for 536 nuclei with A > 40 are shown in Fig. 1201 The dependence of charge radii on the quadrupole 
deformation P has been taken into account for the rare-earth deformed nuclei. 

In the upper left and right panels of Fig. I20( the charge radii for the most stable 159 nuclei with yl > 40 
along the /9-stability line have been analyzed by using the A^^^ and Z^/^ dependence. In the middle left and 
right panels, the same has been done for the data of 536 nuclei with A > 40. Two significant features can be 
observed: (A) On the one hand, the agreement between the data and the calculated results using the Z^^^ 
dependence is much better than that using the A^^^ law, i.e., while there exists a global regular decrease 



38 







50 100 
N 



150 



1.0 1.5 2.0 2.5 3.0 

T1=N/Z 



Figure 21: Left: The rms charge radii versus the neutron number for even-even Ca, Ni, Zr, Sn, Pb isotopes. 
The RCHB calculation with (5-force is represented by open symbols, while the corresponding data is denoted 
by solid symbols. The dashed lines represent the predictions by the A^^^ law with va = 1.228 fm. Right: 
The experimental (solid symbols) and RCHB predicted (dashed lines) coefficient rz = Rc/Z^^^ for the nuclear 
charge radii as a function of isospin quantity r] = N/Z in even-even Ca, Ni, Zr, Sn and Pb isotopes. An 
asymptotic behavior is drawn as a solid line. Taken from Ref. [129]. 

of rA with A, rz nearly remains constant {rz = 1.631(11) fm). The relative rms deviations a for the Z^^^ 
dependence are much less than those for the A^^^ law. (B) On the other hand, though the rms deviation for 
the Z^^^ dependence is significantly reduced, an isospin induced scattering of the data in the middle panels 
in Fig. 1201 can be also observed compared with that in the top panels. In Refs. [38,140,145,146], the isospin 
effects has been considered based on the A^^^ law. As the Z^^'^ dependence can describe the nuclear charge 
radii much better than the A^^^ law, the Z^/^ dependence is a more reasonable starting point for describing 
the isospin dependence of nuclear charge radii. 

To make the isospin dependent Z^^'^ formula for nuclear charge radii be based on more strong foundation 
and be also valid for exotic nuclei, the charge radii of nuclei far from the /3-stability line up to drip lines are 
needed. However, such data is not well available, an alternative is to require that the new isospin dependent 
^1/3 formula should assort with a reliable and microscopic nuclear model such as the fully self-consistent and 
microscopic RCHB theory. 

For even-even Ca, Ni, Zr, Sn and Pb isotopes, the rms charge radii Tc obtained from the RCHB theory (open 
symbols) and the data available (solid symbols) for the even-even Ca, Ni, Zr, Sn and Pb isotopes are given 
in Fig. 1211 The RCHB calculations reproduce the data very well. For a given isotopic chain, an approximate 
linear N dependence of the calculated rms charge radii Tc is clearly seen in Fig. 1211 which shows that the 
variation of r^ for a given isotopic chain deviates from both the simple Z^^^ dependence and the simple A^^^ 
law (denoted by dashed lines in Fig. I21() . Therefore, a strong isospin dependence of nuclear charge radii is 
necessary for nuclei with extreme N/Z ratio. 

In Fig. 1211 the data and RCHB predicted rz = Rc/Z^^'^ for the proton magic isotopes are also presented as 
a function of r/ = N/Z. It is clearly seen that the coefficient rz increases linearly with rj (except some deviations 
due to deformation or shell effects) and the slopes are nearly the same for these isotopic chains. The linear rj 
(or isospin Tz = {N — Z)/2) dependence of rz for an isotopic chain may be understood as the effect of the 
first order perturbation correction of nuclear wave function due to an isospin Tz dependent interaction [147]. 
Based on the analysis of data in the middle and upper panels of Fig. [201 and RCHB prediction in Fig. 1211 the 
following isospin dependent Z^^^ formula for nuclear charge radii has been proposed [129]: 

Rc = rzZ^/'^[l + b{7]-7]*)], v = N/Z, (96) 

where r]* is rj = N/Z for the nuclei along the /3-stability line which can be directly extracted from the nuclear 
mass formula [16], rz = 1.631(11) fm as obtained in upper right panel of Fig. 1201 and b = 0.062(9) obtained 
from the least square fitting. 

The analysis of the available data using Eq. (|96|) with rz and b thus obtained is displayed in the lower right 
panel of Fig. I^UJ The data are reproduced better by Eq. than by Eq. Description of the data with 
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an isospin dependent A^^'^ formula, Rc = 1.228 A^^^ [1 + 0.029(t/ — rj*)] is also shown in the lower left panel of 
Fig- EDI for comparison. The parameter for the isospin term is also obtained by least square fitting of the data. 
It can be seen that this formula describes the data less successfully than Eq. 

In Refs. [140, 145], a modified A^^^ formula with additional ~ {N — Z)/A and ~ 1/^4 terms made some 
successes on describing nuclear charge radius. In addition, this formua was extended with ~ A~'^ terms [38], 
which are considered to account for the isospin dependence of charge radii in light nuclei. Based on the fact 
that RpZ^^^^ is approximately a constant, where Rp is the rms radius of proton density distribution, a formula 
for nuclear charge radius was proposed in Refs. [38,141,142], R^ = ((rpZV3)2 + o.64)i/2. Detailed analysis of 
these expressions show: (i) The Z^^^ dependence is a better approximation to available Rc data than the A^^^ 
one;(ii) The simple Z^^^ dependent formula 1)941) gives nearly the same result as the formula suggested in Refs. 
[141,142]; (iii) The new isospin dependent Z^/^ formula is better than the two-parameter isospin dependent 
^1/3 formula but less better than the three-parameter one. A better agreement of the isospin dependence Z^^^ 
formula of charge radii can be achieved by fitting both rz and b simultaneously. While in Eq. ()96() . the same 
rz is used as that in Eq. (|94)) . which are derived only from the most stable nuclei. Thus the second term in 
Eq. ()96() mainly corresponds to the isospin dependence of charge radii for the nuclei away from the /3 stability 
line. It is expected that the modified Z^^^ formula (|96j) will become more useful with more and more data 
obtained for the nuclei far from the /3-stability line. 

4.2.4 Density distribution and cross section 

The measurement of interaction cross sections plays an important role in exploring the properties of exotic 
nuclei because it can provide information of the nuclear density distribution and nuclear radii. Based on 
the measurement of interaction cross section with radioactive beams at relativistic energy, novel and entirely 
unexpected features appear: e.g., the neutron halo and skin as the rapid increase in the measured interaction 
cross-sections in the neutron-rich light nuclei. Systematic investigation of interaction cross sections for an 
isotope chain or an isotone chain can provide a good opportunity to study the density distributions over a wide 
range of isospin. From the mean field point of view, the nucleon density distributions in an atomic nucleus 
determine the mean potential. Since the single particle levels and the shell structure are directly determined by 
the mean fields, the study of the isospin dependence of shell structure and the mean potential, which become 
highly diffuse near the particle drip line due to the large diffuseness in the density distribution, is also crucial to 
understand properties of unstable nuclei. Particularly the study of the surface diffuseness will provide another 
mean to understand exotic nuclei because the surface diffuseness can be strongly related with the spin-orbital 
splitting which is mainly determined by the diffuseness in the mean field potentials and could be measured 
exp er iment ally. 

The density distribution and shell structure Systematic RCHB calculations have been carried out for 
all the nuclei in Na isotopes with mass number A ranging from 17 to 45 with parameter set NLSH [45]. The 
calculated binding energies Eb and the interaction cross sections with the Glauber Model are shown in Fig. 1221 
The calculated binding energies Eb are in good agreement with the empirical values [89]. The resonance states 
of ^"^Na and ^^Na (with a positive Fermi energy) are exactly reproduced. ^^Na is bound but unstable against the 
proton emission, reproducing the experimental observation. The neutron drip line nucleus has been predicted 
to be ^^Na. The difference between the calculations and the empirical values for the stable isotopes is from the 
deformation [45]. 

The density distributions are examined and the relation between the development of halo and shell effect 
within the RCHB theory is also studied in detail for Na isotopes [45]. For example, the density distributions 
for both the proton and the neutron in Na isotopes are given in Fig. 1231 As seen in the upper part of Fig. 1231 
the change of the neutron density is as follows: for the nucleus with less neutrons (N), the density at the 
center is low and it spreads only to some smaller distance. With the increase of N, the density near the center 
increases due to the occupation of the 2si/2 level, so does the development of neutron radius. In the lower part 
of Fig. 1231 the proton density shows different behavior. The surface is more or less unchanged because of the 
Coulomb Barrier, but with the increase of N, the density of the center decreases due to the slight increasing at 
the tail. This is due to the symmetry energy. But as the density must be multiplied by a factor Anr'^ before 
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Figure 22: Upper: The interaction cross sections cr/ of Na isotopes on a carbon target at 950yl MeV: the open 
circles are the result of RCHB and the available experimental data (A = 20 — 23, 25 — 32) are given by solid 
circles with error bars. The dashed line is a simple extrapolation based on the RCHB calculation for ^^~'^^Na. 
Lower: Binding energies for Na isotopes, the convention is the same as the upper part, but the RCHB result 
for particle unstable isotopes are indicated by triangle. Taken from Ref. [45]. 
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Figure 23: The neutron (upper) and proton (lower) density distributions in Na isotopes. The same figures but 
in logarithm scale are given as inserts to show the tail part of the density distribution. Taken from Ref. [45]. 
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Figure 24: The same as Fig. |221but for matter density distribution. The upper part is given in logarithm scale 
(the radius is labeled at the top of the figure) and the radius vq at which /?(ro) = 10""^ for different isotopes is 
given as inserts. Taken from Ref. [45]. 



the integration in order to give the fixed Z, the big change in the center does not infiuence the outer part of 
the proton distribution very much. 

The matter densities for the even Na isotopes are given in Fig. I24l where the shortest tail in the total density 
occurs for ^^~^^Na, the most stable ones. For either the proton or the neutron rich sides, the tail density 
increases monotonically. The tail (r > 10 fm) of the proton rich nuclei is mainly due to the contribution of the 
proton and that of the neutron drip line nuclei is mainly due to the contribution of the neutron. Compared 
with the neutron-rich isotopes, the proton distribution for nucleus with less N has higher density at the center, 
lower density in the middle (2.5 < r < 4.5fm), a larger tail in the outer (r > 4.5 fm) part. 

It is found that there is interesting connections between the matter distribution and the single particle level 
distribution (see Fig. I25() : after ^^Na, as it is a sub-closure shell for the 1^5/2; the neutrons are filled in the 
2si/2 arid 1^3/2- So the tail of the density for ^^Na is two order of magnitude larger than ^^Na at r = 10 fm, 
while the tail of the density from ■^''Na to '^^Na is very close to each other. But as more neutrons are filled 
in, the added neutrons are filled in the next shell I/7/2) 2^3/2 and 2pi/2- So again two order of magnitude's 
increase has been seen from ^^Na to ^'^Na, and then a gradual increase after ^^Na. So it becomes clear that the 
rapid increase in the cross section is connected with the filling of neutrons in the next shell or sub shell. In the 
inserts in Fig. 1241 the radius rg at which p(ro) = 10~^ fm~^ is given as a function of the mass number to see 
the relation between the shell effect and matter distribution more clearly. It is very interesting to see a slight 
decrease of rg from proton drip line to ^^Na. The tail of ^^Na is the smallest. From ^^Na to '^^Na, one sees a 
almost constant tq. After a jump from '^^Na to ^^Na, a rapid increasing tendency appears again. 

In Figs. 23, the microscopic structure of the single particle energies in the canonical basis [35,42] is given. 
In the left panel of Fig. 1251 the single particle levels in the canonical basis for the isotopes with an even neutron 
number are shown. Going from A = 19 to ^ = 45 we observe a big gap above the N = 8, N = 20 major shell , 
and = 14 sub-shell. The N = 28 shell for stable nuclei fails to appear, as the 2^3/2 and 2pi/2 come so close 
to 1/7/2- When N > 20, the neutrons are filled to the levels in the continuum or weakly bound states in the 
order of 1/7/2, 2^3/2, 2pi/2, and 1/5/2- In the right part, the occupation probabilities in the canonical basis of 
all the neutron levels below E = 10 MeV have been given for ^^Na to show how the levels are filled in nuclei 
near the drip line. 
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Figure 25: Left: Single particle energies for neutrons in the canonical basis as a function of the mass number. 
The dashed line indicates the chemical potential. Right: The occupation probabilities in the canonical basis 
for ^^Na. Taken from Ref. [45]. 



The surface diffuseness and the spin-orbital splitting The RCHB theory has also been used to sys- 
tematically study the tin isotopes from proton drip line to neutron drip line [124]. The spin-orbital splitting 

21 + 1 

versus A for the whole isotope chain is given in Fig.l^Hfor the neutron spin-orbital partners (1^3/2, ^d^/2)^{^97/2^ ^99/2) 
(1^11/2,1^13/2), (l?'i/2,lP3/2), (1/5/2,1/7/2), and (I/19/2, 1/111/2), and the proton spin-orbital partners (1^3/2, 14/2), 
and (1/5/2, 1/7/2)- It is very interesting to see that the spin-orbital splitting for the neutron and proton is very 
close to each other, at least for (1^3/2,1^5/2), and (1/5/2,1/7/2) cases. The splitting decreases monotonically 
from the proton drip line to the neutron drip line for all the partners. To see the underlying reason of this 
behavior, it is very helpful to examine the origin of the spin-orbital splitting in the Dirac equation. The equa- 
tion of motion for the nucleon moving in scalar and vector potentials could be de-coupled and reduced for the 
upper component and the lower component, respectively. If it is reduced in the lower component, it will be 
related with another interesting topic — the pseudo-spin symmetry discussed in next subsection. As for the 
spin-orbital splitting, the Dirac equation can be reduced for the upper component as the following: 

.d^ 1 d{2M-V + S) d ij 

^dr2 E + 2M-V + S dr dr^ 



k{1 + k) 1 Kd{2M -V + S) 



r 



■2 E + 2M - V + S r dr 



where 



{E + 2M -V + S){E-V - S)d^[r), (98) 

I I i = i-i/2. 



The spin-orbital splitting is due to the corresponding spin-orbital potential 

1 K d{2M -V + S) 



(100) 



E + 2M-V + Sr dr 

with some proper normalization factor. It is seen that the spin-orbital splitting in RMF is energy dependent 
and they depends on the derivative of the potential 2M — V + Sas well as the particle distribution. Therefore 

one can introduce the so-called spin-orbital potential: Vis = ^ ; ~ — ^- The Vis for ^^'^Sn, ^^'^Sn, and 

r dr 

^^''Sn, are given in the left panel of Fig. 1271 The Vis for both proton and neutron are almost the same, as the 
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Figure 26: The neutron spin-orbit splitting Ei 
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versus the mass number A in tin 



isotopes for neutron (1^3/2, 1^5/2), (197/2,199/2), (l^ii/2, 1*13/2), (lPi/2, IP3/2), (I/5/2, 1/7/2), and (I/19/2, 1^1/2) 
orbital and proton (1^3/2,1^5/2), and (1/5/2,1/7/2) orbital, respectively. Taken from Ref. [124]. 



potential V — S is a, big quantity (~ 700 MeV), the isospin dependence in the spin-orbital potential could be 
neglected. Therefore the proton and neutron Vis are the same in the present model. That is the reason why the 
spin-orbital splitting for the neutron and proton is very close to each other in Fig. 1261 From ^^'^Sn to ^'^''Sn, the 
amplitude of Vis decreases monotonically due to the surface diffuseness. Neutron potentials V — S for ^'^"Sn, 
^^''Sn/^'^Sn, ^^''Sn, ^■*''Sn,^^'^Sn, ^^^Sn, and ^'''^Sn are given in the right panel of Fig. [27|and its inserted figure 
gives the radii Rq at which V — S = 100 MeV as a function of the mass number. As seen in the above equations, 
the spin-orbital splitting is related with the derivative of the potential V — S . The surface diffuseness happens 
for both the vector and scalar potential, as well as for V — S and V + S. 

For the decline of the spin-orbital splitting, it comes from the diffuseness of the potential or the outwards 
tendency of the potential. The diffuseness of the neutron potentials V + S are given in the left panel of 
Fig. EHl for ii°Sn, ^^Osn, isogn, i^OSn, isogn, leog^^^ no^^_ ^^^^ ^.^at the depth of the potential 

decreases monotonically from the proton drip line to the neutron drip line and the surface of the potential 
moves outwards. The inserted figure gives the radii Rq at which V + S = —10 MeV as a function of the mass 
number. The proton potentials V + S for ^^'^Sn,^^'^Sn, ^^^Sn, ^^'^Sn,^^'^Sn, ^^'^Sn, and ^™Sn are given in the 
right panel of Fig. 1281 With the increase of the neutron number, the proton potentials are pushed towards 
outside as well due to the proton-neutron interaction. 



Interaction cross sections For medium and high energy reactions, the interaction cross section cjint can be 
calculated with the Glauber model which is based on the individual nucleon-nucleon collisions in the overlap 
volume of the high energy colliding nuclei. In the framework of this model, the nucleus-nucleus reaction cross 
section can be obtained from nucleon-nucleon reaction cross section. Consider the collision of a projectile 
nucleus P on a target nucleus T when the nuclei P and T are situated at an impact parameter b relative to 
each other, the total reaction cross section can be written as 

a = 2TT j bdb{l-T{b)), (101) 

where T(b) is the transparency function describing the probability that the projectile will pass through the 
target without interacting. T{b) is given by 



T(6) = exp 



J pp(hp, zp)dhpdzp pT(hT, ZT)dhTdzT t{h — bp -|- bx) 



(102) 



where gnn is the nucleon-nucleon reaction cross section suitably averaged over the interacting n-n, p-p, and 
n-p pairs [148, 149]. The profile function t(b) for the NN scattering is taken as a delta function for zero-range 
nuclear interaction. 
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Figure 27: Left: The derivative of the neutron potentials V{r) — S{r) for ^^'^Sn, ^^^Sn, and ^^'^Sn; Right: The 
neutron potentials V{r) — S{r) for ^'^'^Sn, ^^''Sn/^'^Sn, ^'^'^Sn, ^^'^Sn,^^''Sn, ^^''Sn, and ^'''^Sn. In order to examine 
the surface diffuseness more clearly, the radii Rq at which V{Ro) — S{Ro) = 100 MeV has been given as an 
inserted figure. Taken from Ref. [124]. 




Figure 28: Left: The neutron potentials V{r) + S{r) for ^^^Sn, ii°Sn/20Sn, ^^Ogn, i^Og^^isogj^^ leog^^^ ^^^^ 
^'^''Sn. In order to examine the surface diffuseness more clearly, the radii Rq at which ^(-Ro) + S{Rq) = —10 
MeV has been given as an inserted figure; Right: The proton potentials V{r) + S{r) for ^^'^Sn,^^''Sn, ^^''Sn, 
^^'^Sn,^^'^Sn, ^^'^Sn, and ^^''Sn. In order to examine the surface diffuseness more clearly, the radii Rq at which 
V{Ro) + S{Ro) = — 10 MeV has been given as an inserted figure. Taken from Ref. [124]. 
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Figure 29: The proton density distributions predicted by RCHB for the nuclei C, N, "O and 
i6-25p logarithm scale. Taken from Ref. [46]. 

To compare the cross section directly with experimental measured values, the densities Pn,p{f) of the target 
^^C and the Na isotopes obtained from RCHB (see Fig. I28j) were used in the Glauber model calculation. The 
cross sections for reaction of Na isotopes at 950^ MeV on ^^C have been compared with the experimental 
values [150] in the upper part of Fig. [221 The agreement between the calculated results and measured ones are 
fine. The cross section below ^^Na changes only slightly with the neutron number, which means the proton 
density has played an important role to remedy the contribution of less neutron. From ^^Na to the neutron 
drip line, a gradual increase of the cross section has been observed. After ^^Na, although no data exist yet, a 
relatively fast increase has been predicted. In short, the measured cross section shows similar behavior as that 
of RCHB. 

Charge changing cross section Systematic investigation of interaction cross sections for an isotope chain 
or an isotone chain can provide a good opportunity to study the density distributions over a wide range 
of isospin [45, 150]. However the contribution from proton and neutron are coupled in the measurement of 
interaction cross section. To investigate possible differences in proton and neutron density distributions, a 
combined analysis of the interaction cross section and other experiment on either proton or neutron alone are 
necessary. The charge-changing cross section which is the cross section for all processes which result in a change 
of the atomic number for the projectile can provide good opportunity for this purpose. In Ref. [151], the total 
charge-changing cross section cTcc for the light stable and neutron-rich nuclei at relativistic energy on a carbon 
target were measured. The charge changing cross section dec are calculated by using the Glauber Model and 
the proton density of projectiles and total density of target provided by RCHB theory [46]. 

The proton density distributions predicted by RCHB for the nuclei io-22q^12-24]\j-^14-26q ^^^j^ i6-25p ^^^^ 
given in Fig. EHl in logarithm scale. The change in the density distributions for each isotope chain in Fig. EHl 
occurs only at the tail or in the center part. The density is multiplied by a factor 47rr^ before an integration in 
order to get proton number or radii, therefore the change of density in the central region does not matter very 
much. What is important is the density distribution in the tail part. Compared with the cases of neutron-rich 
isotopes, the proton distribution of nuclei with less has higher density in the center, lower density in the 
middle part (2.5 < r < 4.5fm), a larger tail in the outer part (r > 4.5 fm) which gives rise to the increase of 
and (Tec for the proton rich nuclei as will be seen in following figures. 

With the proton density from RCHB calculations for the projectile used for pp{hp,zp) in Eq. ()102|) . the 
total charge-changing cross sections Ucc of the nuclei ^''"^■^C, ^■^~^^N, 14~26q i6-25p carbon target at 
relativistic energy are calculated and given in Fig. I3U1 The results of the Glauber model are represented by 
open circles and available data [151] by solid ones with error bars. The agreement between the calculation and 
the measurement is very good. 

The charge changing cross sections change only slightly with the neutron number except for proton-rich 
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Figure 30: The total charge-changing cross sections dec of the nuclei io~22q^12-24j\j^14-26q ^^^^ i6-25p ^ 
carbon target at relativistic energy. The open circles are the result of RCHB and the available experimental 
data are given by solid circles with error bars. Taken from Ref. [46]. 



nuclei. This indicates that the proton density plays an important role in determining the charge-changing cross 
sections dec- A gradual increase of the cross section can be observed towards the neutron drip line. However, 
the big error bars of the data cannot help to conclude anything here yet. It is shown clearly that the RCHB 
theory, when combined with the Glauber model, can provide reliable description for not only interaction cross 
section but also charge changing cross section. From comparison between Fig. EO] and Fig. 1171 one finds similar 
trends of variations of proton radii and of charge changing cross sections for each isotope chain which implies 
again the important role that proton plays in determining the charge-changing cross sections. 



4.3 Single particle levels and pseudospin symmetry 

The concept of pseudo-spin is that the single particle orbitals with j = 1+1/2 and j = {1+2) — 1/2 lie very close in 
energy and can therefore be labelled as pseudo-spin doublets with quantum number n = n— 1,/ = / — 1, and s = 
s = 1/2. This concept is originally found in spherical nuclei [28,29] , but later proved to be a good approximation 
in deformed nuclei as well [152-154]. It is shown that pseudo-spin symmetry remains an important physical 
concept even in the case of triaxiality [155]. The origin of pseudo-spin is proved to be connected with the 
special ratio in the strength of the spin-orbit and orbit-orbit interactions [156, 157] and the unitary operator 
performing a transformation from normal spin to pseudo-spin space have been discussed [156-158]. However, 
it is not explained why this special ratio is allowed in nuclei. 

The relation between the pseudo-spin symmetry and the RMF theory was first noted in Ref. [159], in which 
Bahri et al found that the RMF explains approximately the strengthes of spin-orbit and orbit-orbit interactions 
found by non-relativistic calculations. More details have been given in Refs. [153,160], in which it was suggested 
that the origin of pseudo-spin is related to the strength of the scalar and vector potentials. Ginocchio first 
revealed that pseudo-orbital angular momentum is nothing but the "orbital angular momentum" of the lower 
component of the Dirac wave function [30]. 

In the Dirac equation of the nucleon, when the scalar potential ^(r) and the vector potential V{v) are equal 
in amplitudes but opposite in sign, i.e., S{t) + V{t) = 0, or more generally, ^[^(r) -|- V{Y)\/dr = 0, there is an 
exact pseudospin symmetry in single particle spectra [30-32]. Under these conditions, there should be some 
special relations between the four components of the Dirac wave functions of the pseudospin doublets which 
have been used to test the pseudospin symmetry in realistic nuclei [161-164]. 

The pseudo spin symmetry in deformed nuclei has also been extensively studied in the framework of the 
RMF model. The conditions for an exact pseudo spin symmetry and the relationships for the lower components 
of the Dirac eigenfunctions were investigated in Refs. [165-168]. The relationships between the upper and lower 
components of the two states in the pseudo spin doublets have been studied thoroughly for realistic deformed 
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21 + 1 



for Zr and Sn isotopes. From left to right, the pseudo-spin partners correspond 



to (1(^3/2, 2si/2), (1/5/2,2^3/2), (157/2,24/2) and (2^3/2,851/2), respectively. Taken from Ref. [31]. 



relativistic eigenfunctions [164,169]. 

There are also many other investigations on the pseudo-spin symmetry in the framework of RMF the- 
ory [170-177]. A thorough review is certainly beyond the purpose of the present paper. In the following, the 
pseudo-spin symmetry in RCHB theory will be briefly reviewed. 

From the radial Dirac equation for spherical nuclei, similar as Eq. (|98j) . one can derive the radial equation 
for the lower components: 



^dr^ ^ E 



1 d{V + S) d 



V-S)Fl^{r), 



V -S dr 
{E + 2M -V + S){E- 



k(1 



1 



K d{V + S) 



E-V -Sr 



dr 



(103) 



where «;(k — 1) = /(/ -|- 1). It is clear that one can use either Eq. (|103j) or equivalently Eq. (|98)) to get the 
eigenvalues E and the corresponding eigenfunctions. In Section 14.2.41 Eq. ()98|) has been used to discuss the 
spin-orbital splitting in connection with the corresponding spin-orbital potential. If the Eq. (|103j) is used instead 

1 Kd{V + S) 



and the pseudo spin-orbital potential (PSOP) term. 



-, is neglected, then the eigenvalues 



E-V -Sr dr 

E for the same I will degenerate. This is the phenomenon of pseudo-spin symmetry observed in Refs. [28,29]. 
As shown in Ref. [43] , the particle levels for the bound states in canonical basis are the same as those by solving 
the Dirac equation with the scalar and vector potentials from RCHB. Therefore Eqs. H98() and ()103|) remain 
the same in canonical basis even after the pairing interaction has been taken into account. 

In Eq. ()103() . the term which splits the pseudo-spin partners is simply the PSOP. The hidden symmetry 
for the pseudo-spin approximation is revealed as d(y + S)/dr = 0, which is more general and includes V + 
5 = discussed in [30] as a special case. For exotic nuclei with highly diffuse potentials, diV + S)/dr ~ 
may be a good approximation and then the pseudo-spin symmetry will be good. But generally, d[V + 
S) /dr = is not always satisfied in the nuclei and the pseudo-spin symmetry is an approximation. However, if 
. 1 Kd{V + S). k{1-k) 



'E-V -Sr 



dr 



the pseudo-spin approximation will be good. Thus, the comparison of 



the relative magnitude of the pseudo centrifugal barrier (PCB), 



k(1 — k) 



, and the PSOP can provide us with 



some information on the pseudo-spin symmetry. In Ref. [31], with the scalar and vector potentials derived 
from a self-consistent RCHB calculation, the pseudo-spin symmetry and its energy dependence have been 
discussed in realistic nuclei. Afterwards, the pseudo-spin splitting for exotic nuclei in Zr and Sn isotopes has 
been studied [32]. 
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Figure 32: The comparison of the effective pseudo centrifugal barrier (PCB) (E — V — S) 



k{k — 1) 



(dashed 



hnes and dot-dashed hnes) and the effective pseudo-spin orbital potential (PSOP) — 



Kd{V + S) 



dr 



(solid line ) in 



The dashed lines are for 1^3/2 and 2si/2 5 and the 



arbitrary scale for 6^3/2 (upper) and S1/2 (lower) in -"^^^Zr 
dot-dashed lines are for 2^3/2 and 3si/2- The inserted boxes show the same quantities, but the ordinate is 
magnified and the abscissa is reduced to show the behaviors of the effective PCB and the effective PSOP near 
the nuclear surface. Taken from Ref. [31]. 



K- - — K- - 

The redefined pseudo-spin orbital splitting AE = — —z — — versus the average single particle 

+ (I + l)-£'i,-=/'_i/2 

energy E = — for the bound pseudo-spin partners in Sn and Zr isotopes are plotted 

in Fig. 1311 In both isotopes, a monotonous decreasing behavior with the decreasing binding energy is clearly 
seen. The pseudo-spin splitting for 35^/2 and 2^3/2 is more than 10 times smaller than that of the 2sx/2 and 
1^3/2- As far as the isospin dependence of the pseudo-spin orbital splitting is concerned, the splitting in Sn 
isotopes gives a monotonous decreasing behavior with the increasing isospin. Particularly for 2sii2 and 1^3/2 
partners, the pseudo-spin splitting in ^^''Sn is only half of that in ^^Sn. Just as one expects [31], the pseudo-spin 
symmetry in neutron-rich nuclei is better. In Zr isotopes, although the situation is more complicated (e.g., 
the effect of the deformation which is neglected here), the pattern is more or less the same, i.e., a monotonous 
decreasing behavior with the decreasing binding energy and a monotonous decreasing behavior with the isospin. 
From these studies, the pseudo-spin symmetry remains a good approximation for both stable and exotic nuclei. 
A better pseudo-spin symmetry can be expected for the orbital near the threshold, particularly for nuclei near 
the particle drip line. 

To understand why the energy splitting of the pseudo-spin partner changes with different binding energies 
and why the pseudo-spin approximation is good in RMF, the PSOP and PCB should be examined carefully. 
Unfortunately, it is very hard to compare them clearly, as the PSOP has a singularity aX E ^ {V + S). As one 

is only interested in the relative magnitude of the PCB and the PSOP, the effective PCB, {E — V — S)— — 2 — ^> 

and the effective PSOP, ^ — 7-^- — ^ are introduced for comparison. They correspond to the PCB and the 

r dr 

PSOP multiplied by a common factor E — V — S respectively. 

The effective PSOP does not depend on the binding energy of the single particle level, but depends on the 
angular momentum and parity. On the other hand the effective PCB depends on the energy. Comparing these 
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Figure 33: The comparison of the effective pseudo centrifugal barrier (PCB) [E — V — S) ^ (dashed 

hues) and the effective pseudo-spin orbital potential (PSOP) ^ — — ^ (solid line ) multiplied by the square 

r dr 

of the wave function F of the lower components in arbitrary scale for (upper) and Si/2 (lower) in ^^'^Zr. 
Taken from Ref. [31]. 



two effective potentials one could see the energy dependence of the pseudo-spin symmetry. They are given in 
Fig. EH for si/2 (lower) and (upper) of ^^°Zr in arbitrary scale. 

The pseudo-spin approximation is much better for the less bound pseudo-spin partners, because the effective 
PCB is smaller for the more deeply bound states. This is in agreement with the results shown in Fig. 1311 The 
effective PSOP and the effective PCB are also given as inserts in Fig. Ellin order to show their behavior near 
the nuclear surface. 

The effective PCB (dashed lines or dot-dashed lines) and the effective PSOP (solid lines) multiplied by the 
squares of the lower component wave function F{r) are given in Fig.inHl for 2si/2 (upper left), 3si/2 (lower left), 
1^3/2 (upper right), and 2^3/2 (lower right) of ^^''Zr in arbitrary scale. The pseudo-spin approximation is much 
better for the less bound pseudo-spin partners, because the effective PCB is smaller for the more deeply bound 
states. This is in agreement with the results shown above. It is clear that the contribution of the effective PCB 
(dashed lines or dot-dashed lines) is much bigger than that of the effective PSOP (solid lines). Generally the 
effective PSOP is two orders of magnitude smaller than the effective PCB. The energy difference between the 
orbital j = I + 1/2 and the orbital j = I — 1/2 is always negative with some exceptions. Further investigations 

reveal that the integration of '^'^ |Fp over r gives the splitting of the pseudo-spin partners, whose sign 

dr 

will decide the normal splitting or the reverse. 



4.4 Halos and giant halos 

Since the experimental discovery of neutron halo phenomena in ^^Li [13], the study of exotic nuclei has become 
a very challenging topic in nuclear physics. Experimentally more and more halo nuclei have been observed [6- 
8, 11, 12, 178-181]. On the theoretical side very different models have been used to describe these phenomena. 
Three-body calculations based on an inert core of ^Li and two outside neutrons treat the translational invariance 
properly, but they neglect polarization effects. Full shell model calculations contain all the configuration mixing, 
but they are limited to relatively small configuration spaces in the oscillator model. In the mean field description, 
the large radius observed experimentally in ^^Li has been interpreted by the fact that filling in more and more 
neutrons in the nuclear well, the Fermi surface for the neutrons comes close to the continuum limit. Due to 
the small single- neutron separation energy, the tails of the two wave functions of the last filled orbital (lpi/2) 
reach very far outside of the nuclear well and a neutron halo is formed [13]. Although this simple interpretation 
is based on the mean field picture, several microscopic theoretical investigations within self-consistent mean 
field models [182-185] have failed. In fact, it would be a strange accident if the last occupied neutron level 
in ^^Li, the lpi/2 orbit, would be so close to the continuum limit that the tail of the last two uncorrelated 
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neutron wave functions reaches so far out as would be necessary to reproduce the large experimental radius 
observed. Therefore, in the non-relativistic scheme Bertsch et al [182] and Sagawa [184] introduced an artificial 
modification to the potential in order to reproduce the small separation energy in their mean field calculations 
using Skyrme interactions. In this way the authors were able to reproduce qualitatively the observed trends, 
although some discrepancies still remained. 

Koepf et al [183] were the first to investigate the neutron halo in the RMF theory for both spherical and 
axially deformed cases. They found large deformations for the lighter Li isotopes, but a spherical shape for ^^Li 
and as in the non-relativistic investigations the binding energy of the lpi/2 too large so as to reproduce a 
neutron halo. Zhu et al [185] improved the result of the RMF calculations by applying a similar modification 
to the potential as in Ref. [184] in order to adjust the proper size of the halo. 

The first self-consistent microscopic description of the halo in ^^Li was achieved by using the RCHB theory 
where the pairing correlations are taken into account by a density dependent force of zero range and the 
halo was reproduced in a self-consistent way by using the scattering of Cooper pairs to the 2si/2 level in the 
continuum [36]. 

In Ref. [121], pairing correlations and the coupling to particle continuum states are described by finite range 
two-body forces and the finite element method are used in the coordinate space discretization of the coupled 
system of Dirac-Hartree-Bogoliubov integro-differential eigenvalue equations, and Klein-Gordon equations for 
the meson fields [71]. The authors of Ref. [121] performed calculations for the isotopic chains of Ne and C 
nuclei and found evidence for the occurrence of neutron halo in heavier Ne isotopes. The RHB theory is also 
applied to the description of the properties of light nuclei in isotopic chains of C, N, O, F, Ne, Na, and Mg 
with large neutron excess in Ref. [186]. The RHB model predicts the formation of neutron skin and eventually 
a neutron halo in Ne and Na due to the quasi-degeneracy of the triplet of states 1^1/2, 2p3/2) and 2pi/2 which 
are close to the Fermi surface in neutron rich nuclei. 

To describe medium-heavy nuclei far from the /3 stability line, the main problem of nuclear structure 
models is the extrapolation of effective forces to nuclei with extreme isospin values because most of these nuclei 
are still not accessible experimentally and therefore one could not compare theoretical results with empirical 
data. In order to make predictions for medium-heavy nuclei at the neutron drip line, one has to test available 
effective interactions in detailed calculations of the properties of neutron-rich nuclei for which a comparison 
with experimental data is possible. This has been done in Ref. [187] where the RHB theory is applied in the 
description of the ground-state properties of Ni and Sn isotopes. In a comparison with available experimental 
data, it is shown that the NL3 -|- Gogny DIS effective interaction provides an excellent description of binding 
energies, neutron separation energies, and proton and neutron rms radii. The results indicate that this choice 
of model parameters might also be valid for nuclei with more extreme isospin values, i.e., medium- heavy nuclei 
at the drip lines. 

We note that the RHB model solved in a discretized coordinate space is available only for spherical systems 
up to now. The deformed RHB theory solved in a harmonic oscillator basis [188] is not applicable for halo 
nuclei due to the localization property of the harmonic oscillator potential. One should work in two or three 
dimensional spatial lattice or instead, one can improve the asymptotic behavior of the HO wave function [64,65] 
or adopt other basis which has a correct asymptotic behavior [73]. 

In the following, we will briefly review the RCHB description of halo phenomena and the prediction of giant 
halo in light and medium-heavy nuclei. 

4.4.1 Neutron halo in light nuclei 

The ground state properties of Li isotopes with mass numbers A = 6 to A = 11 were investigated by using the 
RCHB theory with the effective interaction NL2 [36]. Pairing is neglected for the three protons and the strength 
Vq of the pairing force H84|) for the neutrons is determined by a calculation in the nucleus ^Li adjusting the 

corresponding pairing energy -— Tr Ak to that of a RHB-calculation in an oscillator basis using the finite range 

part of the Gogny force DIS of Ref. [119] in the pairing channel. Good agreement with experimental values is 
found for the total binding energies and the radii of the isotope chain ^Li to ^^Li. The matter radius shows a 
considerable increase when going from the nucleus ^Li to ^^Li. In contrast to the earlier mean field calculations 
of Refs. [182, 184, 185] these results are obtained without any artificial modifications of the potential. 
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Figure 34: Calculated and experimental density distribution in ^^Li and ^Li. The solid line shows the result 
of ^^Li while the dashed line corresponds to the calculation of ^Li. The shaded area gives the experimental 
results with error bars. Taken from Ref. [36]. 
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Figure 35: The mean field potential S + V for protons (l.h.s.) and neutrons (r.h.s.). The chemical potential is 
given by a dashed line. The energy levels in the canonical basis are indicated by horizontal lines with various 
lengths proportional to the occupation of the corresponding orbit. Taken from Ref. [36]. 

In Fig. I^the corresponding density distribution for the neutrons in the isotopes ^Li and ^^Li is shown. 
It is clearly seen that the increase of the matter radius is caused by a large neutron halo in the nucleus ^^Li. 
Its density distribution is in very good agreement with the experimental density of this isotope shown with its 
error bars by the shaded area. 

In order to understand the microscopic structure of this halo, in Fig. E^lthe mean field S{r) + V{r) for the 
protons and neutrons together with the energy levels in the canonical basis is shown. The Fermi level for the 
neutrons is very close to the continuum limit in close vicinity to the and to the z^2si/2 level. The length 

of these energy levels in Fig. [33 is proportional to the corresponding occupation. One clearly finds that pairing 
correlations cause a partial occupation of both the i'lpi/2 and the z^2si/2 level, i.e. a scattering of Cooper pairs 
to the continuum. 

In contrast to many previous investigations, in the RCHB model, the halo is not formed by two neutrons 
occupying the lpi/2 level very close to the continuum limit, but is formed by Cooper-pairs scattered mainly 
in the two levels lpi/2 and 2si/2- This is made possible by the fact that the 2sx/2 comes down close to the 
Fermi level in this nucleus and by the density dependent pairing interaction coupling the levels below the Fermi 
surface to the continuum. The RCHB theory provides a much more general mechanism for halo occurring in 
nuclei. One needs only several single particle levels with small orbital angular momenta and correspondingly 
small centrifugal barriers close to, but not necessarily directly at, the continuum limit. 
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Figure 36: Canonical basis single-particle neutron levels as functions of the number of neutrons. The spectrum 
is calculated for Ne isotopes. Taken from Ref. [121]. 

In Ref. [121], the isotopic chains of Ne and C nuclei were studied with the RHB model with finite range 
pairing interaction in the coordinate space. The formation of neutron skin and eventually a neutron halo in Ne 
was found to be mainly due to the quasi-degeneracy of the triplet of states lf7/2, 2p3/2, and 2pi/2- In Fig- OHl 
the neutron single-particle spectrum in the canonical basis is displayed for Ne isotopes. The energies of levels 
in the continuum decrease with increasing neutron number (cf. Fig. The shell structure dramatically 

changes at > 22. The triplet of states lf7/2, 2p3/2, and 2pi/2 approaches zero energy, and a gap is formed 
between these states and all other states in the continuum. The Fermi level uniformly increases toward zero for 
A'^ = 22. Between N = 22 and = 32, the Fermi level is practically constant and very close to the continuum. 
The addition of neutrons in this region of the drip does not increase the binding. Only the spatial extension of 
neutron distribution displays an increase. At = 32 the Fermi energy becomes slightly positive, and heavier 
isotopes are not bound any more. The finite range pairing interaction promotes neutrons from the lf7/2 orbital 
to the 2p levels. Since these levels are so close in energy, the total binding energy does not change significantly. 
Because of their small centrifugal barrier, the 2p3/2 and 2pi/2 orbitals form the halo. 

4.4.2 Giant halos in Zr isotopes 

In all of the halos observed so far, one has only a very small number of nucleons, namely one or two outside 
of the normal core. In order to study the influence of correlations and many-body effects it would be very 
interesting to find also nuclei with a larger number of neutron distributed in the halo. Following the successful 
description of the halo nucleus ^^Li, the RCHB theory was used to study and predict giant neutron halos for 
Zr-isotopes close to the neutron drip line [43]. It is formed by up to six neutrons outside of the ^^^Zr core with 
the magic neutron number N = 82. 

From the upper panel of Fig. |37| which presents the rms radii of the protons and neutrons for the Zirconium 
isotopes, one finds a kink for the neutron rms radius at the magic neutron number A^ = 82. This kink can 
be understood more clearly by considering the microscopic structure of the underlying wave functions and the 
single particle energies in the canonical basis which is shown in the lower panel of Fig. 1371 Going from A^ = 70 
to A^ = 100, there is a big gap above the I/111/2 orbit. For A^ > 82, the neutrons are filled to the levels in the 
continuum or weakly bound states in the order of 3p3/2, 2/7/21 '^Pi/2j 2/5/2 and 1/19/2- 

The neutron chemical potential is given in this Figure by a dashed line. It approaches rapidly the continuum 
already shortly after the magic neutron number N = 82 and it crosses the continuum at A^ = 100 for the nucleus 
to ^-^^Zr. In this region the chemical potential is very small and almost parallel to the continuum limit. This 
means that the additional neutrons are added with a very small, nearly vanishing binding energy at the edge of 
the continuum. The total binding energies E for the isotopes above ^■^^Zr are therefore almost identical. This 
has been recognized already in Ref. [189] in RMF calculations using the BCS approximation and an expansion 
in an oscillator basis, which is definitely not reliable for chemical potentials so close to the continuum limit. 

In the upper panel of Fig. |3Slthe density distributions for neutrons and protons in the nucleus ^^''Zr are 
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Figure 37: Upper: rms radii for neutrons and protons in Zr isotopes close to the neutron drip line as a function 
of the mass number A. Lower: single particle energies for neutrons in the canonical basis. The dashed line 
indicates the chemical potential. Taken from Ref. [43]. 
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Figure 38: Upper: neutron and proton density distribution in ^ Zr. Dashed lines indicate calculations for 
different values of the box size R and the dashed-dotted line give the neutron distribution for the core ^^^Zr. 
Lower: relative contributions of the different orbits to the full neutron density as a function of the radius. The 
shaded area indicates the total neutron density in arbitrary units. Taken from Ref. [43]. 
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Figure 39: The occupation probabilities in the canonical basis for various Zr isotopes with mass number A as 
a function of the single particle energy. The chemical potential is indicated by a vertical line. For A > 122 the 
number Nf^ of neutrons in the halo is also shown. Taken from Ref. [43]. 

shown. A very large box size is needed to describe the halo properly. For R =30 fm the neutron density is 
reliably reproduced only up to r = 25 fm, where it has decreased to 10~^ fm~^. The full line in the upper 
panel of Fig. EHl is an asymptotic extension to infinite box size. On the other hand, the density distribution 
inside the nucleus is reproduced properly even for small values of the box size. The relative contributions pnij 
of the different orbits characterized by the quantum numbers nlj with respect to the total neutron density pn 
are displayed in the lower panel of Fig. 1381 The shaded area gives the total neutron density in arbitrary units. 
The halo is formed essentially by contributions from three orbits 3p3/2, 3pi/2) and 2/7/2- The most inner part 
of the halo (7 < r < 9 fm) the 2/7/2 orbit plays the dominant role. As can be seen in the lower part of Fig. 1381 
this orbit is slightly below the chemical potential and to the continuum limit in this nucleus. Further outside 
(10 < r < 15 fm) its relative contribution is strongly reduced because of the larger centrifugal barrier felt by 
the / = 3 orbit. In this region the orbit 3p3/2, which has nearly the same position as the 2/7/2 orbit, takes 
over. Because of the smaller orbital angular momentum, this orbit feels a reduced centrifugal barrier. For even 
larger distances from the center (r > 15 fm) its relative contribution is somewhat reduced and the 3pi/2 orbit 
gains importance. The 3p3/2 and the 3pi/2 levels feel the same centrifugal barrier, but the latter is situated 
directly at the continuum limit and therefore it is more loosely bound than the other two orbits. 

The occupation probabilities in the canonical basis of all the neutron levels near the Fermi surface, i.e. 
in the interval —20 < E < W MeV are presented in Fig. |3H for all the Zr-isotopes between A = 108 and 
A = 140. The chemical potential is indicated by a vertical line. For the mass numbers A < 122 the chemical 
potential lies several MeV below the continuum limit {E = 0) and there is only very little occupation in the 
continuum {E > 0). The nucleus ^^^Zr has a magic neutron number and no pairing. As the neutron number 
goes beyond this closed core, the occupation of the continuum becomes more and more important. Adding up 
the occupation probabilities v"^ for the levels with E > one finds in the continuum 2 particles for N = 84, 
4 for = 86, 6 iV = 88, roughly 3 for A^ = 90, roughly 4 for A^ = 92, roughly 3 for A^ = 94, roughly 4 for 
A^ = 96, roughly 5 for A^ = 98, and roughly 6 for A^ = 100 (where the neutron drip line is reached). 

Considering the very extended neutron rms radii of these systems and estimating the number of valence 
neutrons which would fit into the same volume outside the ^■^^Zr core (shown in Fig. |2Hlby a dashed dotted 
line) if they would be packed with normal neutron density, one finds for ^^^Zr 24 neutrons instead of 12 and for 
^^'^Zr 34 neutrons instead of 18. This means that the density outside of the ^■^^Zr core is reduced by roughly a 
factor of 2. This phenomenon is therefore clearly a neutron halo and not a neutron skin. As one sees from Fig. 
3 the tail of the density is proportional to exp(— 2\/4m5r) with S ^ 0.4 MeV. This fact is also emphasized by 
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the extremely low 2n-separation energies (< 0.5 MeV) of these systems. This phenomenon is called a Giant 
Halo because of the large number of particles in the halo region [43].. 

4.4.3 Giant halos in lighter nuclei 

Although the giant halos were predicted in Zr isotopes, it is very difficult for experimentalists to confirm them 
because the exotic Zr isotopes with > 82 are too heavy to be synthesized by the RIB facilities at present. 
Thus it is valuable to investigate the giant halo phenomena in lighter nuclei which are experimentally accessible 
today. Ground state properties of all the even-even O, Ca, Ni, Zr, Sn, Pb isotopes ranging from the proton 
drip line to the neutron drip line were investigated [44, 126]. The two neutron sepation energies are shown in 
Fig. I14l where one finds the S2n values for exotic Ca isotopes near neutron drip line are very close to zero in the 
large mass region, i.e., 52„ w 2.06, 1.24, 0.76, 0.43, 0.21, 0.04 MeV for A = 62, 64, 66, 68, 70, 72 isotopes, 
respectively. If one regards ^''Ca as a core, then the several valence neutrons occupying levels above the N = 40 
sub-shell for these A > 60 exotic nuclei are all weakly bound and can be scattered easily to the continuum 
levels due to the pairing interaction, especially for ^^~'^^Ca. This case is very similar to S2n in Zr isotopes 
with N > 82 (see Fig. I14() [43] . Note that for Sn isotopes in the vicinity of the drip line [124], however, S2n 
decrease rather fast with the mass number and for Ni isotopes [42] , S2n are quite large (~ 2MeV) and undergo 
a sudden drop to negative value at the drip line nucleus. Such behavior of S2n in Ca isotopes with A > 40 
indicates that there exist giant halos in these nuclei, just as what happens in Zr isotopes with A > 82 [43]. 
From Fig. 1181 one also finds similarity between Ca and Zr isotopes: the neutron radius increases suddenly with 
neutron number increasing in the Ca chain above A = 40 and in Zr chain above A^ = 82. Together with the 
analysis of nucleon density distributions, single particle levels and the occupation probabilities in the canonical 
basis, it is predicted that giant halo may exist in Ca isotopes with A>60 [44,126]. 

Compared with giant halos in Zr chain, the halos in Ca isotopes will be much easier to access with today's 
radioactive facilities. This means that the giant halos in the Ca isotopic chain are challenging but promising 
for experimentalists. In order to make the conclusion more solid and exclude the effects from possible error in 
S2n (~ 1 MeV), in Fig. 1401 the two- neutron separation energies S2n in drip line region for even- neutron Ne, Na, 
Mg, and Al chains (left panel) and for Ar, K, Ca, Sc and Ti chains (right panel) are presented. Open symbols 
represent the values calculated from the RCHB theory with the NLSH parameter set while corresponding filled 
symbols represent the data available. It should be noted that the calculations are performed by keeping the 
shape spherical. Nevertheless it can be seen from the figure that all these S2n curves are almost parallel to 
each other and the zero S2n line. Although the calculated S2n value may differ from the data by 1 or 2 MeV, 
there must be one or more isotope chains in which many nuclei have the S2n values very close to zero thus are 
good candidates for giant halos. 

Experimentally much efforts have been made to search the drip line for heavier elements. New experiments 
have shown that both ^'^Na and ^^Ne are bound while ^^Ne and ^^Na are missing [190, 191], as expected from 
the predictions in Ref. [70]. In Fig. I4()| ^"^Na and ^^Ne lie in the drip line region and approach the predicted 
giant halo area. Therefore much more experimental effort should be devoted to measuring their masses and 
extending the isotope chains for the possible giant halos. 

4.5 Halos in hypernuclei 

The production mechanisms, spectroscopy, and decay modes of hypernuclear states have been the subject 
of many theoretical studies. The relativistic Hartree Bogoliubov model in coordinate space with finite range 
pairing interaction [71,121] has also been applied to describe A hypernuclei with a large neutron excess [192]. It 
is found that the inclusion of the A hyperon does not produce excessive changes in bulk properties but shifts the 
neutron drip by stabilizing an otherwise unbound core nucleus at drip line. For example, the nucleus ^^Ne was 
unbound without the A (see Fig. ^2 (a))- The presence of the strange baryon stabilizes the otherwise unbound 
core and ^^^'^Ne becomes bound (see Fig. (b)). The microscopic mechanism through which additional 
neutrons are bound to the core originates from the increase in magnitude of the spin-orbit term in presence of 
the A particle [192]. 

Motivated by experimental knowledge of A-N interaction and theoretical understanding on giant halo [43, 
44, 126], it is interesting to observe the possible appearances of halos in hyper exotic nuclei as what have been 
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Figure 40: Left: The two-neutron separation energies S2n for even-neutron Ne, Na, Mg, and Al nuclei in 
the drip hne region. Open symbols represent the values calculated with the RCHB theory with the NLSH 
parameter set, while corresponding solid symbols represent the data available. The horizontal line at 2MeV 
denotes the upper limit for the possible occurrence of halos. Right: for even Ar, K, Ca, Sc, and Ti nuclei in 
the drip line region. Taken from Ref. [126]. 
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Figure 41: lf-2p single-particle neutron levels in the canonical basis for the Ne (a), and Ne-|-A (b) isotopes. 
The dotted line denotes the Fermi level. Taken from Ref. [192]. 



done in Ref. [47]. In Fig. 1421 two neutron separation energies S2n for ordinary nuclei, single-A hyper nuclei 
and double-A hyper nuclei of Ca isotopes labelled by A = 0, A = 1 and A = 2, respectively, from the proton 
to neutron drip line are presented. It is clear that one or two A hyperons lower the neutron Fermi level a little 
bit but keep the neutron shell structure unchanged. Therefore, the neutron drip line is pushed outside from 

= 52 in ordinary isotope chain to A = 54 in hyper isotope chain. Meanwhile, giant halo due to pairing 
correlation and the contribution from the continuum still exist in Ca hyper nuclei similar as that in ordinary 
Ca isotopes [44, 126]. This is a slight but rewarding step for exploring the limit of drip line nuclei on the basis 
of the giant halo. For more details, see Ref. [47]. 

Apart from neutron halos in hyper nuclei, as A hyperon is less bound than the corresponding nucleon in 
nuclei, it is worth investigating the existence of hyperon halos. In Ref. [194], a lambda halo was tentatively 
suggested in a RMF plus BCS model. However, the normal BCS method will have unphysical solution by 
involving baryon gas [42], such prediction needs further study. Therefore hyper carbon isotopes are studied 
by RCHB theory in Ref. [193]. In Fig. 1431 the baryon density distributions in logarithmic scales versus the 
radius r in hyper nuclei j^^C, and 3^C are presented, where the solid, dot-dashed, long-dashed and dotted 
lines represent density distributions for the lambda, neutron, proton, and corresponding core, respectively. It 
can be seen that up to two A hyperons added to the core ^^C, the nucleon density distributions remain the 
same and hyperon density distributions at the tail are comparable with those of the nucleons. An intriguing 
phenomenon appears in g^C where the hyperon density distribution has a long tail extended far outside of its 
core (denoted by dotted line), as shown in the right panel of Fig. 1431 This is a signature of hyperon halo. 
It is due to the weakly bound state Ip^/o in which has a small hyperon separation energy and a density 
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Figure 42: The two-neutron separation energies S2n in even-N Ca isotopes versus the neutron number N. The 
upper panel is for ordinary nuclei, the middle one for single-A hyper nuclei, and the lower one for double-A 
hyper nuclei. Taken from Ref. [47]. 
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Figure 43: Density distributions for A (solid), neutron (dot-dashed), proton (long-dashed) and corresponding 
core (dotted) in J^^C, and g^C in logarithmic scales. Taken from Ref. [193]. 

distribution with long tail. 

5 Magic numbers for superheavy nuclei 

One of the fundamental and persistent questions in nuclear science is the exploration for the limits of charge and 
mass that a nucleus can attain and the creation of nuclei with masses and charges much larger than those we 
are familiar with, i.e., superheavy nuclei. The properties of superheavy nuclei are important to understand not 
only the nuclear structure, but also the structure of the stars and the evolution of the universe. However, after 
lots of efforts for both higher luminosities and better detection efficiencies as well as the impressive progresses 
on the synthesis of superheavy nuclei, the borders of the upper-right corner in the nuclear chart still remain 
unknown. 

Due to fission resulting from the balance of the attractive nuclear surface tension and the repulsive Coulomb 
force, the superheavy nuclei cannot exist in the classical charged liquid-drop model. The stability of superheavy 
nuclei is mainly determined by shell effects. At the magic proton or neutron numbers 2, 8, 20, 28, 50, and 
82, as well as N=126 for neutrons, nuclei have higher stability and abundance compared with their neighbors. 
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Particularly, the highest stability is observed in the case of the doubly magic nuclei. Amongst other special 
properties, the doubly magic nuclei are spherical. A semblance of the same would work also for superheavy 
nuclei, if there were magic numbers in this region. Consequently, these nuclei will be guarded against a faster 
decay by fission and have more opportunities to be bound. Therefore, it is important to find out the regions 
in the (Z, N) plane where the shell effects are strong and the long lifetimes of the superheavy nuclei can be 
expected. 

The pioneering works on the superheavy elements are a series of macro-microscopic nuclear shell model cal- 
culations performed over more than 30 years to determine nuclear masses and energy surfaces, which suggested 
the existence of a group of relatively stable nuclei separated in neutron and proton numbers from the known 
heavy elements by a higher instability when approaching the closed spherical shells Z=114 and N=184 [195-198]. 
This group of nuclei become to be known as the island of superheavy elements and the verification of their 
existence is one of the most challenging topics in world-wide heavy ion research facility. 

Experimentally, however, although the superheavy elements up to Z=116 are synthesized or claimed to be 
synthesized [199-212] by different heavy ion reaction types including the cold, warm and hot fusion reactions, 
they are not examples of the originally thought island of superheavy elements yet. 

There is no consensus among different theories with regard to the center of the island of superheavy nuclei. 
Based upon phenomenological models such as finite-range droplet model (FRDM), the shell closures were 
predicted at Z=114 and A^=184 [213]. Additionally, the FRDM also predicts larger shell gaps at Z=104, 
106, 108, 110 and at Ar=162, 164 [213]. In Nilsson-Strutinsky scheme, a similar pattern of deformed nuclei 
have been predicted about Z=108 and iV=162 as in FRDM [214,215]. However, the main obstacle is the 
question whether the macroscopic approaches which apply to the region of /^-stability line can be extrapolated 
to the superheavy nuclei. Recently, microscopic calculations [87, 216-226] are also attempted in describing 
the superheavy nuclei. In the framework of RHB theory, calculations with a finite-range pairing force of 
Gogny interaction Dl and effective interaction NLSH show that Z=11A and iV=160, iV=166, A''=184 exhibit 
stability compared to their neighbors and indications for a doubly magic character at Z=106 and iV=160 are 
also observed [217]. The Skyrme Hartree-Fock (SHF) method with interactions SkP and SLy7 predicts magic 
numbers at .Z=126 and N=184, and also predicts the increased stability due to the deformed shell effects at 
N = 162 [218]. Considering non-relativistic SHF effective interactions SkM*, SkP, SLy6, Skll, SkIS, SkI4 
and relativistic mean field (RMF) effective interactions PL-40, NLSH, NL-Z, TMl, the doubly magic spherical 
nuclei i84ll4, 172I2O and 184I26 based on two-nucleon gaps 52p and 52n are given [219]. The uncertainty on 
the magic numbers lies in the uncertain strength of the spin-orbit coupling in non-relativistic approaches or 
the effective interactions in relativistic models. The prediction of magic numbers in superheavy nuclei remains 
a challenge for the nuclear models, which have to be rigorously tested by a wide variety of nuclear properties 
throughout the periodic table. 

As the RCHB formalism allows for the proper description of the coupling between the bound states and the 
continuum by the pairing force and has shown a remarkable success in the description of nuclei with unusual 
N/Z ratio [36,42-44,46], it is suitable to apply the RCHB theory for superheavy nuclei. In principle, only a 
calculation in a large multidimensional deformation space can definitively decide the appropriate ground-state 
shape. However, as the traditional superheavy nuclei are expected to be spherical and are located on the 
nuclear chart around a spherical doubly magic nucleus next to f^i-fb, the RCHB theory with the assumption 
of spherical shape can be applied for the search. Therefore here we will mainly focus on the magic numbers 
in spherical superheavy nuclei. Some crucial observables in searching shell closures and magic numbers are 
taken into account, including two-nucleon separation energies, two-nucleon gaps, the shell correction energies, 
as well as the pairing energies and the effective pairing gaps. At the end of this section, a brief discussion on 
the stability of the doubly spherical nuclei against deformation and the main progress on the investigation of 
the deformed superheavy nuclei will be given, as well as the a— decay half-lives of superheavy nuclei. 

5.1 Two nucleon separation energies 

Traditionally, the magic nuclei are characterized by higher stability than their neighboring ones and the magic 
numbers can be revealed by the differences of the binding energy, i.e. the one- or two-nucleon separation 
energies [16]. Therefore, the nucleon separation energies are good starting points to investigate the magic 
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numbers in superheavy region. Due to the absence of odd-even effects, the two-nucleon separation energy 5*2^ 
(52n) is better to quantify shell effects than the single-nucleon separation energy Sp (5„). For an isotonic 
(isotopic) chain, the S2p {S2n) decrease with the proton (neutron) number and its sudden jump indicates the 
occurrence of a proton (neutron) magic number. 
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Figure 44: The two-proton(left panel) and two-neutron (right panel) separation energies, S2p and S2n, as a 
function of mass number A obtained by RCHB calculation with PKl effective interaction. Taken from Ref. [48]. 

In Ref. [48], the two-proton and two-neutron separation energies have been investigated in the superheavy 
region with proton number Z ranging from 100 to 140 and neutron number N from {Z + 30) to {2Z + 32) 
by the RCHB theory with NLl, NL3, NLSH, TMl, TW99, DD-MEl, PKl and PKIR. As an example, the 
results with PKl are shown as a function of mass number A in Fig. 1441 For the proton in left panel, the large 
gaps appear at Z = 106, 120, 132, 138, while the usual Z=11A does not appear. Furthermore, all the gaps 
are strongly neutron dependent. On the neutron side in right panel, the relatively larger gaps can be seen at 
A^=138, 172, 184, 198, 228, 238, 252, 258, and 274, and also the size and the shape of the gaps differ notably. 

From similar investigations with other effective interactions [48], it is found that the magic proton numbers 
Z=120, 132, and 138 are common while Z=106 is observed only for NL3, NLSH, TW99, PKl, and PKIR. 
Furthermore the gaps do not exist for all isotopes. On the neutron side, S2n show distinguishable gaps at 
A^=138, 172, 184, 198, 228, 238, 252, 258, and 274 with all the effective interactions. Apart from these, there 
are also other gaps which appear only for some effective interactions, e.g., N = 164 for NLSH, TW99, DD-MEl, 
PKl and PKIR, and N = 216 for NLSH, TW99, PKl, and PKIR. For the magic numbers marked either by 
the common gaps or the interaction-dependent gaps, the shell quenching phenomena, i.e., the gaps at A^=184, 
198 (NLl), 216 (NLSH, TW99, PKl, and PKIR), 228 (NLl, NL3, NLSH, TMl, PKl, and PKIR), 238, 252, 
and 258 (NLSH, TMl, PKl, and PKIR) appear only for certain Z, are observed. 



5.2 Two-nucleon shell gaps 

The changes of the two-nucleon separation energies can be quantified by the second difference of the binding 
energies, i.e., the two-nucleon shell gaps: 

52p{N,Z) = S2p{N,Z)-S2p{N,Z + 2) = 2Eb{N,Z)-Eb{N,Z + 2)-Eb{N,Z-2) (104) 
52n{N,Z) = S2n{N,Z)-S2n{N + 2,Z) = 2EB{N,Z)-EBiN + 2,Z)-EB{N-2,Z). 

A pronounced peak of the two-nucleon shell gaps corresponds to a drastic change of the two-nucleon separation 
energies and indicates the shell closure [227]. So far, the two-nucleon shell gaps have been extensively used to 
be an indicator for the magic number [219,221] and to analyze the shell quenching phenomenon [228]. 

The two-proton gaps 52p from the RCHB calculations for even-even nuclei with Z=102 - 138 as a function 
of Z are shown in Fig. I45L The two-proton gaps S2p with the same are connected as a curve. A peak at 
certain Z in the curve suggests the existence of magic proton number. The sharpness of the peaks represent the 
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Figure 45: The two-proton shell gaps, 62p, as a function of proton number obtained by RCHB calculation with 
effective interactions NLl, NL3, NLSH, TMl, TW-99, DD-MEl, PKl, and PKIR, respectively. Taken from 
Ref. [48]. 
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Figure 46: The two-neutron shell gaps, as a function of neutron number obtained by RCHB calculation 
with effective interactions NLl, NL3, NLSH, TMl, TW-99, DD-MEl, PKl, and PKIR, respectively. Taken 
from Ref. [48]. 
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goodness of the magic numbers while the quenching effects are associated with the bundle of the curves at the 
certain Z. The size of the gaps of two-proton separation energies S^-p in Fig. correspond to the magnitude 
of the peaks of two-proton gaps in Fig. 1451 It is observed that common magic proton numbers Z=120, 132, 
and 138 exist for all the effective interactions while Z=106 is observed only for NL3, NLSH, TW99, DD-MEl, 
PKl, and PKIR. Furthermore, the peak at Z=114 for NLSH and TW99 and the peak at Z=126 for NLl are 
also observed, though they are not so obvious as that at Z=120. 

Similar figure for two-neutron gaps 62n are shown in Fig. 1461 There are interaction-independent peaks of 
two-neutron gaps 62n at A^=138, 172, 184, 198, 228, 238, 258, and 274 which are consistent with the conclusion 
from the two-neutron separation energies S2n ■ Moreover, small peaks are also observed for A^=154 (NLSH 
and TW99), 7V=164 (NL3, NLSH, TW99, DD-MEl, PKl, and PKIR), and iV=216 (NLSH, TW99, PKl, and 
PKIR). 

5.3 Shell correction energies 

The shell correction energy, representing the overall behavior of the single-particle spectra, is another suitable 
observable to identify the shell closure. This quantity, which is derived from the single-particle spectra by 
the Strutinsky procedure [229], is defined as the difference between the total single-particle energy E and the 
smooth single-particle energy E: 

N(Z) j 

Esheii = E-E=^ei-2 / egie)de, (105) 
where N{Z) is the particle number, Cj is the single-particle energy, A is the smoothed Fermi level determined 

A 

by the particle number equation N{Z) = 2 J g{e)de, and g{e) is the smoothed level density 
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with the smoothing range 7 , the folding function f{x) = —^e ^ P{x), and P{x) the associated Laguerre 



polynomial lI^'^{x'^). The shell correction energy provides an indicator about the deviation in the level structure 
of nuclei away from uniformly distributed ones. A large negative shell correction energy corresponds to shell 
closure at certain nucleon number. 

With the proton single-particle spectra in canonical basis from the RCHB calculations, the shell correction 
energies E^^^^^ for even-even nuclei with Z ranging from 100 to 140 are shown as a function of Z in Fig. 1471 The 
nuclei with the same N are connected as a curve. A deep valley at certain Z on a curve hints the shell closure. 
The valleys at Z=120, 132, and 138 are conspicuous for all the effective interactions. Possible shell closure 
at Z=106 for NLl, NL3, NLSH, TW99, DD-MEl, PKl, and PKIR, at Z=1U for NLSH, TW99, PKl, and 
PKIR, and at Z=126 for NLl are also observed. The shell closure from the shell correction energies for proton 
^ shell consistent with the peaks of the two-proton gaps S2p, though the magnitude of magicity slightly 
differs. Similar to the two-proton gaps 62p, the spread of the valleys of the shell correction energies for proton 
^ shell indicates the quenching of the magic number. In view of the shell correction energies for proton -E'^/^g;;, 
the quenching phenomena of magic proton numbers in RCHB calculation are universal. These facts suggest 
that the magic proton numbers depend also on neutron number N. 

The corresponding shell correction energies for neutron -E^^g;; for even-even nuclei with A^=130 - 312 are 
demonstrated as a function of N in Fig. 1481 As the neutron numbers extend from A^=130 to A^=290, the 
valleys for shell correction energies for neutron E'^^^^ in Fig. ^HJare not so obvious as those for shell correction 
energies for proton E-^^^^^ in Fig. 1471 However, similar conclusions as the two-neutron separation energies S2n 
and the two-neutron gaps 62n for magic neutron numbers can also be drawn. Similar calculation has also been 
done in Ref. [230] with SHE and RMF. 
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Figure 47: The shell correction energies for proton, E^^^^, as a function of proton number obtained by RCHB 
calculation with effective interactions NLl, NL3, NLSH, TMl, TW-99, DD-MEl, PKl and PKIR, respectively. 
Taken from Ref. [48]. 
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Figure 48: The shell correction energies for neutron, -E^^g;;, as a function of neutron number obtained by RCHB 
calculation with effective interactions NLl, NL3, NLSH, TMl, TW-99, DD-MEl, PKl and PKIR, respectively. 
Taken from Ref. [48]. 
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The first valley locates in A^=184 for all the effective interactions and involves N=172 (NL3, NLSH, TMl, 
TW99, DD-MEl, PKl, and PKIR), even iV=164 (NLSH and TW99). And the other valley is the mixture 
of the valleys at N=228 (TW99 and DD-MEl), iV=252 (NLSH, TMl, PKl, and PKIR) and A^=258. These 
fine structures suggest the smearing of magic neutron numbers. Based on the shell correction energies, the 
shell closures are smeared compared with those from the two-nucleon separation energies 5*2^ and S2n and the 
two-nucleon gaps 62p and 62n- However, it is found in Ref. [48] that the magic numbers at Z=120, 132, and 
138 are common but strongly quenched, the shell closures near A^=172 and A^=184, and N=228, N=252, and 
A^=258 are blurred. The other magic numbers from the shell correction energies, such as Z=W6, 114, and 126, 
appear only for some effective interactions. 

In the usual Hartree or HF mean field level, as the nonphysical particle continuum are involved, much 
efforts have been made to avoid the divergence of single particle energy density around the threshold [231-234]. 
Based on the Green's-function approach to the level density, a method of calculating shell corrections has 
been adopted within the SHF and RMF theories to a large-scale survey of spherical shell energies throughout 
the whole landscape of conceivable superheavy nuclei [230]. From SHF with SkP, SLy6, SkI3 and SkI4 and 
RMF with NL3 and NL-Z2, the total shell correction energies Egheii for spherical superheavy nuclei around 
the expected island of stability around Z = 120 and N = 180 are given in Fig. 1491 Instead of narrow stripes 
of large E shell localized around magic numbers for normal nuclei, a wide area of shell stabilization is found 
to spreads over all shell closures predicted. In other words, the predicted magic shells for superheavy nuclei 
are much mellowed, which are consistent with the results from RCHB theory. The disappearance of a familiar 
pattern of magic numbers and the appearance of broad valleys of shell stability are due to the rather large 
single-particle level density and the appearance of many low j shells around the Fermi level [230]. 
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Figure 49: Total shell correction energies, calculated from SHF with SkP, SLy6, SkI3 and SkI4 and RMF with 
NL3 and NL-Z2, for spherical even-even superheavy nuclei around the expected island of stability around Z = 
120 and = 180. The thick solid lines denote two-particle drip lines. Black squares mark nuclei calculated 
to be stable with respect to (3 decay. White color indicates nuclei with positive shell corrections, black color 
denotes nuclei with Egheii beyond -12 MeV. Taken from Ref. [235]. 



5.4 Magic numbers and deformations 

5.4.1 Magic numbers for spherical superheavy nuclei 

After the detailed analysis of the two-nucleon separation energies, S2p and S2n-, the two-nucleon shell gaps, 
52p and 52n-, the shell correction energies, E^^^^^ and E^f^^^, as well as the pairing effects, the magic num- 
bers for superheavy nuclei in the RCHB theory are found to be [48]: Z = 120,132,138 for proton and 
N = 172, 184, 198, 228, 238, 258 for neutron. As many theoretical methods have been adopted to probe the 
superheavy magic numbers and there is no consensus among them with regard to the center of the island of 
superheavy nuclei. It is interesting to briefly review the up-to-date studies of the superheavy magic numbers 
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including the present one: 

• Z=106: As we have seen, this shell closure is strongly supported by the effective forces NL3, NLSH, 
TW-99, PKl, and PKIR, partially by NLl and DD-MEl, but not at all by TMl. We also noted that the 
shell closure at Z=106 depends sensitively on the neutron number and is visible only around N ~ 154 
and ~ 218, while the later can only be clearly seen in the results calculated with TW-99 and NLSH. The 
RHB calculations with the effective force NLSH and the Gogny interaction Dl for the pairing correlation 
predicts a doubly magic character at Z=W6 and A^=160 [217]. In addition, the shell closures at Z=104, 
106, 108, and 110 and A^=162, and 164 which correspond to prolate deformations are predicted by the 
FRDM [213] and the deformed shell closure occurs at Z=108 in Nilsson-Strutinsky scheme [214,215]. 

• Z=114: The shell closures are predicted at Z=114 and A^=184 in the phenomenological shell models [197, 
198,213,214]. However, as we have noted, in the RCHB theory the magicity of Z=\1A is not at all obvious, 
in agreement with other RMF investigations [22,217,219,236,237]. By adopting a new effective interaction 
NL-RAl, the nucleus ig4ll4 is predicted to be the next spherical doubly magic nucleus in a RMF plus 
BCS calculation [238]. Within the SHF approach, only SkI4 prefers the doubly magic nucleus at Z=114 
and A^=184, which is not supported by other Skryme forces such as SkM*, SkP, Skll, SkI3, SLy6, and 
SLy7 [219,221]. 

• Z=120: The shell closure at Z=120 is strongly supported in the RCHB calculations. The doubly magic 
characteristic of X72I2O has been addressed within the RMF theory [219, 221, 234, 238] as well as the 
effective field theory [237] . Most Skyrme HF investigations also suggest that this nucleus is doubly magic 
except for SkP [219,221]. Due to the strong shell quenching effect in the Z=120 isotopic chain, other 
doubly magic candidates with A^=184, 198, 228, 238, 252, 258, and 274 are interaction-dependent and 
need to be investigated in more detail. 

• Z=126: Several Skyrme HF investigations with the effective interactions such as SkP, SkM*, SLy6 and 
SLy7 predict doubly magic system at Z=126 and A^=184 [218,219,221]. However, this conclusion is not 
supported by the RMF calculations although A^=184 might be a neutron magic number. 

• Z=132 and 138: Little attention has been paid to the possible proton shell closures at Z > 130 due 
to their large Coulomb repulsive potentials. However, the shell closures at Z=132 with N ~ 210-240, 
and at Z=138 with N ~ 180-200 and N ~ 258 appear in most relativistic and non relativistic mean field 
models [48,219], . 

Normally the magic numbers correspond to the large energy gaps between the single particle levels and the 
doubly magic nuclei are those with protons and neutrons filled up to the gap. It is therefore interesting to 
study the single-particle level structure in these doubly magic nuclei, or vice versa. 




Figure 50: Spherical single particle spectra in 172120 for protons (left) and neutrons in various mean- field 
calculations as indicated. Taken from Ref. [221]. 

As examples, the spherical single particle spectra in ^^-^120 for protons (left) and neutrons in various mean- 
field calculations are given in Fig. [201 [221]. The occurrence of the shell closure at Z = 120 depends on the 
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amplitude of the spin-orbit splitting of the 3p states above the Fermi level and the 2f levels below the Fermi 
energy. It has also been pointed out that the shell closures may be coupled to the density profiles (and vice 
versa) based on the systematic analysis of density distribution of superheavy nuclei. For ^^^120, a significant 
central depression or central semibubble is obtained in the spherical case [221,239]. Recently, the influence of 
the central depression in the density distribution of spherical superheavy nuclei on the shell structure is studied 
in some detail within the RMF theory [240]. It is found that a large depression leads to the shell gaps at the 
proton Z = 120 and neutron N = 172 numbers, whereas a flatter density distribution favors = 184 and leads 
to the appearance of a Z = 126 shell gap and to the decrease of the size of the Z = 120 shell gap. On the other 
side, the deformation effect has been found to be significant in the calculation of the density distribution in 
this nucleus, namely only a weak central depression in the deformed case for ^^^120 compared to the predicted 
semibubble at the spherical shape [240,241]. 



5.4.2 Stability of the doubly magic superheavy nuclei against deformation 

The magic numbers discussed above are mainly based on the assumption of spherical geometrical configuration. 
However, such assumption is not always true for superheavy nuclei [213-215,242]. It is worthwhile to investigate 
the potential energy surfaces to see the validity of spherical configuration. 

In Fig. EH the potential energy surfaces for '^^'^120, ^°H20, ^18120, ^^^120, 358120, and ^78120 (solid lines) 
obtained from constrained RMF calculations with the pairing correlations treated by BCS approximations 
are displayed. Although these nuclei have the same proton number, their potential energy surfaces are quite 
different from each other as shown in Fig. 1511 For nuclei ^^^120, ^''^120 and ^"^^120, there is an obvious local 
minimum with the spherical configuration (^2 ~ 0, while another local minimum with large deformation /?2 ~ 0.6 
can be also clearly seen. For ^^^120, the spherical local minimum is very shallow, and for ^48120 and ^58120, 
the spherical minimum is hardly seen and a local minimum with P2 ~ 0.25 appears instead. There are also 
local minimum around P2 ~ —0.4 with oblate deformation. In addition to the local minima discussed above, 
we would like to make a few remarks on the absolute minima for these nuclei. For ^^^120 and '^'^^120, the two 
minima with different deformations have almost the same energies, i.e. the so-called "shape coexistence" may 
exist. In particular, the spherical minimum is indeed the absolute minimum for 292 120 {Eb = -2064.3 MeV 
with /?2 ~ vs. Eb = -2063.6 MeV with P2 ~ 0.6), while the ground state of ^''8120 is quite delicate with Eb = 
-2398.7 MeV for /?2 ~ and Eb = -2398.7 MeV for /?2 ~ 0.6. For 3°^120, the absolute minimum at f32 ~ 0.6 
is much deeper (~ 5.5 MeV) than the spherical configuration. In addition for "^"^^120, the absolute minimum 
lies at P2 ~ 0.2, and it is a well-deformed nucleus according to the present calculation. From the RCHB and 
the constrained RMF+BCS calculations, it is shown that spherical doubly magic nuclei may exist, e.g., ^^^120 
and 378120 [48]. 

The role of shell effects in superheavy nuclei can be also seen in Fig. 1511 in which the corresponding 
macroscopic energy, defined as the difference between the binding energy and the total shell correction energy, 
Eb — -Esheih is plotted with dashed lines. It can be seen that without the shell effects, the superheavy nuclei 
hardly exist. As noted long time ago, the shell effects play an essential role to stabilize the superheavy nuclei 
against the fission. 

Apart from axial symmetric case, the influence of other deformation such as triaxiality and reflection 
asymmetry deformations has been investigated as well, e.g., see Ref. [23,226]. Figure EH presents an example 
of potential energy surface of the typical actinide nucleus ^^''Pu, which has often served as a benchmark for 
mean- field models, obtained with the Skyrme force SLy6. A typical doubled- humped fission barrier is well 
shown. The inner barrier explores triaxial degrees of freedom, which reduce the axial barrier by about 3 MeV, 
while the outer barrier explores refiection asymmetric shape [226]. The similar behavior is held for heavier 
nuclear system [220] . It is also noted that though the further deformation degrees of freedom reduce the fission 
barriers obviously, they hardly change the energies of mimina. Therefore, the constraint RMF model with axial 
symmetry remains valid to confirm the spherical configuration of a superheavy nucleus. 
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Figure 51: The total energy, Eb, and the equivalent macroseopic energy, Eb — Es^cW, of 292,304,318,348,358,378 -|^20 
calculated in the constrained RMF theory with effective interaction NL3. Taken from Ref. [48]. 




Figure 52: Double-humped fission barrier of the typical actinide nucleus ^^'^Pu. The dotted line denotes 
an axial and reflection-symmetric calculation, the fuh line denotes a triaxial (inner barrier) and axial and 
reflection-asymmetric calculation (outer barrier). The various shapes along the axial paths are indicated by 
the contours of the total density at vq = 0.07 fm~^. Taken from Ref. [226] 
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5.5 Alpha decay and the half-hves 



The preferred decay mode of shell stabilized superheavy nuclei is a decay. A key quantity for a decay, the Qa 
value is defined as 

Q„(iV, Z) = E{N, Z) - E{N - 2, Z - 2) - E{2, 2), (107) 

which is the energy of the a-particles emitted by radioactive heavy and superheavy nuclei. It is interesting to 
compare the Qa data of the discovered superheavy nuclei with predictions from mean-field models. Such works 
have been done in various of mean field framework, e.g. Refs. [23,222-225,243]. In the recent experiments 
designed to synthesize the element 115 in the ^^^Am+^^Ca reaction at Dubna in Russia, three similar decay 
chains consisting of five consecutive a-decays, and another different decay chain of four consecutive a-decays 
are observed, and the decay properties of these synthesized nuclei are claimed to be consistent with consecutive 
a-decays originating from the parent isotopes of the new element 115, ^^^115 and ^^'^115, respectively [212]. The 
deformed RMF-I-BCS calculation with a density-independent delta-function interaction in the pairing channel 
and the effective interaction TMA is made to study these newly synthesized superheavy nuclei ^^^115, ^^^115, 
and their a-decay daughter nuclei [225]. The calculated a-decay energies and half-lives agree well with the 
experimental values and with those of the macroscopic-microscopic FRDM-I-FY and YPE-I-WS models. 

For an area of enhanced stability, the a-decay half-lives are expected to be longer than its neighbors. The 
half- lives of a-decay can be obtained with Qa value by phenomenological Viola and Seaborg systematics [244] : 

logioTc, = (1.66175Z - 8.5166)Q^^/2 - (0.20228Z + 33.9069) (108) 

where Z is the proton number of the parent nucleus, Tq, in second and the parameters taken from Ref. [245]. 
It should be noted that Eq. ()108() is based on the WKB approximation, and provides only a rather crude 
estimation of since it disregards many structure effects such as deformation, and configuration changes, etc. 
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Figure 53: The a-decay half-lives Tq as a function of neutron number obtained by RCHB calculation with PKl 
effective interaction. Taken from Ref. [48]. 

In Fig. 1531 the half-lives from the RCHB calculation in logarithm scale are plotted as a function of 
neutron number N with the effective interaction PKl. The half-lives corresponding to 1 ns, 1 fis, 1 ms, 1 s, 
1 h, 1 y, and 1 ky are marked by the dashed lines. Each curve in this figure corresponds to an isotopic chain 
in the region with Z ranging from 102 to 140. In Fig. 1531 the half-lives logioT^ increase with A^. The jumps 
of the curves correspond to magic proton numbers Z=W6, 120, 132, and 138, which depend on the effective 
interactions. The peaks of each curve correspond to magic neutron numbers, i.e., A^=172, 184, 228, 238 and 
258. It can be seen that the magic numbers suggested by two- neutron separation energies (two- neutron gaps 
^2n) or two-proton separation energies S2p (two-proton gaps 52p) can also be found here. It should be noted 
that the half-lives Tq, shown here are obtained from the spherical RCHB calcualtions. In fact, both the pairing 



68 



correlation and deformation are essential to reproduce the experiment Qa values [222,224,246]. Calculations 
along this line are necessary. 

6 Summary and perspectives 

This article has been concerned with the relativistic description of the exotic nuclei particularly with the 
Relativistic Continuum Hartree-Bogoliubov (RCHB) theory and its applications. The last two decades have 
witnessed the development of the relativistic mean field (RMF) theory and its successful application in nuclear 
physics. At the very beginning, the question on the clear signal of relativistic effects in atomic nuclei has puzzled 
us for quite a long time. Now after lots of efforts, it become evident that there are at least the following points 
which do favor the relativistic description of nuclei: 

• One may say that the kinetic energy of a nucleon moving in an nucleus is much smaller than its mass. 
However, as far as the mass is concerned, we have to remember that the mass of a nucleus is smaller 
than the sum of the mass of its consistent nucleons, which is naturally obtained in a relativistic model. 
Furthermore there are always a large scalar and a large vector potentials which are both comparable to 
the mass of the nucleons in nucleus. 

• Relativistically the spin-orbit potential comes out naturally. Furthermore the long existing problems on 
the origin of pseudo-spin symmetry were solved. A very well developed spin symmetry in single anti- 
neutron and single anti-proton spectra is found. This spin symmetry in anti-particle spectra and the 
pseudo-spin symmetry in particle spectra have the same origin. 

• In the framework of the relativistic approach, the nucleons interact via the exchanges of mesons and 
photons. From this point of view, the relativistic approach is more microscopic and fundamental in the 
sense of describing the nuclear system at the meson level. 

• In practice, it turns out that the RMF theory can reproduce better the nuclear saturation properties (the 
Coester line) in nuclear matter and the isotopic shifts in the Pb-region, and is more reliable for nuclei far 
away from the line of /3-stability. 

Altogether the relativistic mean field theory is a reliable and effective model for nuclear structure. Therefore, 
with the fast progress in producing the exotic nuclei, it is natural to describe the exotic nuclei based on the 
RMF theory. 

The formalism, the numerical solutions and the effective interactions for the RMF theory as well as its 
application for nuclear matter and neutron stars are presented in Section 2. The numerical solutions of the 
Dirac equation for spherical nuclei, particularly on the Woods-Saxon basis, are discussed in detail. The new 
effective interactions with nonlinear self-coupling meson fields and density dependent meson-nucleon couplings, 
and their influences on properties of nuclear matter and neutron stars are also discussed. 

Pairing correlation plays important roles in open shell nuclei. The conventional BCS method can be 
combined easily with the RMF theory, but it is not justified for exotic nuclei because it could not include 
properly the contribution of continuum states. The formalism and related discussions of the RCHB theory, in 
which the coupling between the bound states and continuum could be included self-consistently, are given in 
Section 3. Discussion on the continuum states, the resonant BCS method and several methods for obtaining 
single particle resonant states are briefly reviewed. 

The study of exotic nuclei has attracted world wide attention due to their large N/Z ratios (isospin) and 
interesting properties such as halo and skin. In Section 4 we reviewed the application of the RCHB theory to 
exotic nuclei, the binding energies, particle separation energies, the radii and cross sections, the single particle 
levels, shell structure, the restoration of the pseudo-spin symmetry, the halo and giant halo, and halos in hyper 
nuclei, etc. For examples, (1) the Li isotopes were investigated and the halo phenomenon in ^^Li was found 
to be the result of the scattering of particle pairs into continua; (2) the giant halos were predicted in Zr, Ca 
and even lighter nuclei when neutron drip line is approached based on the analysis of two neutron separation 
energies, single particle energy levels, the orbital occupation, the contribution of continuum and nucleon density 
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distribution. These predictions are supported to some extent by the recent experimental information; (3) exotic 
phenomena in hyper-nuclei were also studied; etc. 

The RCHB theory has been applied to study superheavy nuclei which is reviewed in Section 5. The possible 
doubly magic superheavy nuclei have been searched within the RCHB theory using a variant of the effective 
interactions. Base on the detailed analysis on the two-miclcon separation energies and two-nucleon 
shell gaps 52p and 52n-, the shell correction energies Egj^gn and Egii^n, as well as the pairing energies -Ep^ir 
^pairi effective pairing gaps Ap and A„, the proton and neutron shell closures have been predicted by using 
the RCHB theory with the effective interactions NLl, NL3, NLSH, TMl, TW-99, DD-MEl, PKl and PKIR. 
Proton numbers Z=120, 132, and 138 and neutron numbers A^=172, 184, 198, 228, 238, and 258 are supported 
by all effective interactions to be magic. The deformation-constrained RMF calculations have been done for 
these doubly magic nuclei candidates and discussions on the deformation of these nuclei are also presented. 

Among many perspectives of the application of the RMF and the RCHB theories, we mention briefly the 
following: 

1. Attempts have been made by extending the present RCHB theory for the axially deformed nuclei. Al- 
though the coupled Dirac equations and Hartrce-Bogoliubov equations with fixed pairing potentials have 
been solved in coordinate space, the simultaneous self-consistent solution of axially deformed RCHB equa- 
tions is proved to be very difficult [247] . Therefore the Woods-Saxon basis has been suggested to solving 
the RMF theory in order to generalize it to study exotic nuclei by solving either the Schrodinger equation 
or the Dirac equation [73]. It combines the advantages of the efficiencies in harmonic oscillator basis and 
the appropriateness for exotic nuclei in coordinate space. For spherical case, the Woods-Saxon basis has 
been proved to be very successful [73]. The development of deformed RCHB in the Woods-Saxon basis 
is in progress. 

2. The Shell-model-Like APproach (SLAP) method has been suggested to treat pairing correlation in RMF 
theory [248] . With the single particle bound states from the RMF calculation and resonant states obtained 
from RMF theory combined with the Analytic Continuation of Coupling Constant (ACCC) method [110, 
111] or the phase shift scattering method [41,117], the SLAP for the pairing can also be used to describe 
exotic nuclei as well. 

3. Since the pion does not contribute in the Hartree level, we did not mention its contribution throughout 
the paper. However, if the Fock term is considered, one must include pion in the Lagrangian. This would 
bring out some complication in solving the equations as well as new and interesting physics. 

4. Many symmetries have been broken in the mean field theory. Restoration of some symmetries could 
be achieved by making projections or corrections. The microscopic center of mass correction for the 
broken translation invariance in the RMF model has been discussed [51]. The single particle states of a 
deformed mean field do not have good total angular momentum, so is the total intrinsic wave function 
of the deformed nucleus. Thus one must project out a state with good total angular momentum from 
the intrinsic wave function. This can be done by the full angular momentum projection method [35] or 
approximately in a simple way [249]. Recent works on the restoration of rotation symmetry [250,251] 
are very interesting. Possible extensions and applications for chiral doublet bands and nuclear isomeric 
states are expected [252-258] . 

5. In most of the RMF calculations, the no-sea approximation is applied. Therefore the anti-nucleon degrees 
of freedom are not considered. It would be quite interesting to investigate these degrees of freedom and 
the polarization of the vacuum as well as abnormal nuclei containing anti-baryons [33,47,193,259,260]. 
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